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We introduce a method, dubbed the flux-fusion anomaly test, to detect certain anomalous symme¬ 
try fractionalization patterns in two-dimensional symmetry enriched topological (SET) phases. We 
focus on bosonic systems with Z 2 topological order, and symmetry group of the form G = U(l) x G', 
where G' is an arbitrary group that may include spatial symmetries and/or time reversal. The 
anomalous fractionalization patterns we identify cannot occur in strictly d = 2 systems, but can 
occur at surfaces of d = 3 symmetry protected topological (SPT) phases. This observation leads 
to examples of d — 3 bosonic topological crystalline insulators (TCIs) that, to our knowledge, have 
not previously been identified. In some cases, these d = 3 bosonic TCIs can have an anomalous 
superfluid at the surface, which is characterized by non-trivial projective transformations of the 
superfluid vortices under symmetry. The basic idea of our anomaly test is to introduce fluxes of the 
U(l) symmetry, and to show that some fractionalization patterns cannot be extended to a consis¬ 
tent action of G' symmetry on the fluxes. For some anomalies, this can be described in terms of 
dimensional reduction to d = 1 SPT phases. We apply our method to several different symmetry 
groups with non-trivial anomalies, including G = U(l) x iff and G = U(l) x if , where if and if 
are time-reversal and d = 2 reflection symmetry, respectively. 


I. INTRODUCTION 

Following the theoretical prediction [1-6] and experi¬ 
mental discovery [7, 8] of time-reversal invariant topolog¬ 
ical band insulators, it has become clear that symmetry 
plays a rich and varied role in topological phases of mat¬ 
ter. New families of symmetric topological phases have 
been identified theoretically, and significant strides have 
been made in the classification and characterization of 
such phases. Much of the recent progress, with some im¬ 
portant exceptions, has focused on systems with internal 
(or, on-site) symmetry, such as time reversal, U(l) charge 
symmetry and SO (3) spin symmetry. For example, free- 
fermion topological insulators and superconductors with 
internal symmetry have been fully classified [9, 10]. Sub¬ 
sequent work identified the symmetry protected topolog¬ 
ical (SPT) phases, some of which are strongly interacting 
generalizations of topological insulators that do not ad¬ 
mit a free-electron description [11-16]. 

Less attention has been paid to the role of crystalline 
space group symmetry in topological phases, especially 
in the setting of strongly interacting systems. Of course, 
such symmetry is common and varied in real solids, in 
contrast to a relatively small number of realistic inter¬ 
nal symmetries. Therefore, with an eye toward eventual 
experimental realizations of new topological phases, it is 
important to develop theories of such phases with crys¬ 
talline symmetry [17-29]. To accomplish this task, new 
theoretical approaches are needed, as some of the exist¬ 
ing tools to classify and characterize topological phases 
are limited to internal symmetry. 

In this paper, we consider two-dimensional (d = 2) 
topologically ordered systems, where crystalline and 
other symmetries play a non-trivial role via their action 


on anyon quasiparticle excitations [17, 19, 30-33]. Such 
systems are said to be in symmetry-enriched topological 
(SET) phases. We introduce a method, the flux-fusion 
anomaly test , which allows us to show that some puta¬ 
tive SET phases cannot exist in strictly two dimensions. 
However, such states can exist as surfaces of d = 3 SPT 
phases. Our method allows us to identify new exam¬ 
ples of d = 3 SPT phases dubbed bosonic topological 
crystalline insulators (TCIs), which are outside the scope 
of existing theoretical approaches, via their surface SET 
phases. Bosonic TCIs in d = 3, named after electronic 
TCIs [29], are SPT phases where the protecting symme¬ 
try includes both U(l) and the space group symmetry of 
a clean d = 2 surface. These states are interesting not 
only in the context of spin or boson systems, but as a 
possible stepping stone toward understanding electronic 
TCIs with strong interactions, and we hope our results 
can spur more progress in this direction. 

The jumping off point for our approach is a consid¬ 
eration of symmetry fractionalization patterns in d = 2. 
Provided we assume symmetry does not permute anyon 
species, the action of symmetry fractionalizes into an ac¬ 
tion on individual anyons, hence the term symmetry frac¬ 
tionalization. The classic example is the fractional charge 
of Laughlin quasiparticles in fractional quantum Hall liq¬ 
uids [34]. We refer to a complete description of symme¬ 
try fractionalization in a topologically ordered system as 
a symmetry fractionalization pattern. Distinct patterns 
of symmetry fractionalization - including for crystalline 
symmetry - have been classified [17, 19, 31], and the sym¬ 
metry fractionalization pattern is a universal property of 
a SET phase [19]. 

A symmetry fractionalization pattern may be anoma¬ 
lous, which means that it cannot occur in a strictly d = 2 
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Symmetry 

G = U(l) x G' 

Vison fractionalization 
classes \H 2 (G',Z 2 )} 

Anomaly-negative vison 
fractionalization classes (J\T) 

d = 3 SPT phases distinguished 
by anomaly test (<S) 

Anomaly 

type 

U(l) x Zl 

Z 2 

Zi 

Z 2 

1 

U(l) x Ll 
(reflection) 

Z 2 

Zi 

Z 2 

2 

U(l) x pm 
(translation & 
parallel reflection) 

(Z 2 ) 4 

(Z 2 ) 2 

(Z 2 ) 2 

2 

(U(l) x Zl) x pi 
(translation only) 

(Z 2 ) 4 

(Z 2 ) 2 

(Z 2) 2 

3 

(U(l) x Zl) x pm 

(Z 2 ) 8 

(Z 2 ) 3 

(Z 2 ) 5 

2,3 

U(l) x pAmm 
(square lattice) 

(Z 2 ) 6 

(Z 2 ) 3 

(Z 2 ) 3 

2 

(U(l) xZ l) xpAmm 

(z 2 ) 10 

(Z 2 ) 4 

(Z 2 ) 6 

2,3 


TABLE I. Summary of results. Each row is a distinct symmetry group, given in the first column. The last column indicates 
the type or types of anomalies that appear, as described in the text. The meaning of the other columns is discussed in the 
text. Zi denotes the trivial group. In all these cases we consider Z 2 gauge theory whose gauge charge e carries half U(l) charge 
while the gauge flux m carries zero charge. 


system, but is instead realized at the surface of a d = 3 
symmetry-protected topological (SPT) phase [35-38]. In 
this case, we say we have a surface SET phase. SPT 
phases [11-16] have an energy gap, lack spontaneous sym¬ 
metry breaking, and, upon weakly breaking whatever 
symmetries are present, are in the trivial phase; that is, 
the ground state wave function can be adiabatically con¬ 
tinued to a product state when symmetry is explicitly 
broken. It follows that SPT phases lack bulk excitations 
with non-trivial braiding statistics. Instead, edge or sur¬ 
face properties are generally non-trivial; for d = 3 SPT 
phases, one possibility is to have a surface SET phase 
with anomalous symmetry fractionalization. 

While a number of results have been obtained on 
anomalous symmetry fractionalization of internal sym¬ 
metry [35-38], generalization to incorporate crystalline 
symmetry is not straightforward. Our approach, the flux- 
fusion anomaly test, is a method to test for anomalous 
symmetry fractionalization for symmetries of the form 
G = U(l)xG', where G' is an arbitrary group that may 
include crystalline symmetry. We focus on bosonic sys¬ 
tems, such as spin models or systems of bosons. We note 
that some results on anomalous reflection symmetry frac¬ 
tionalization have recently appeared in Ref. 28. We also 
note that the “monopole tunneling” approach developed 
in [36] and used in [37] is closely related but not equiva¬ 
lent to the flux-fusion anomaly test, as discussed further 
in Sec. VIII. 

The basic idea of the flux-fusion anomaly test is to start 
with a symmetry fractionalization pattern for a d = 2 
SET phase, to introduce fluxes of the U(l) symmetry, 
and then to determine whether the fractionalization pat¬ 
tern can be extended to an action of G' symmetry on 
the U(l) fluxes. Sometimes this is impossible, signaling 
anomalous symmetry fractionalization. These considera¬ 
tions only depend on the fusion rules of fluxes and anyon 
excitations, hence the name for the anomaly test. We 
emphasize that we do not need to consider flux thread¬ 


ing or flux insertion as a dynamical process. 

We implement this idea by gauging a subgroup Z„ C 
U(l), and studying the resulting theory. Gauging sym¬ 
metry has been employed to study SPT phases, where 
different phases can be distinguished using the statistics 
of excitations in the gauged theory [39]. Here, the gauged 
theory is itself a SET phase with G' symmetry. We are 
able to show that some symmetry fractionalization pat¬ 
terns are anomalous by studying the action of G' sym¬ 
metry on the anyons of the gauged SET phase. 

We primarily consider symmetries of the form G = 
U(l) x G S p ace and G = (U(l) x Zl) x G space , where 
Z 2 is time reversal and G space is a d = 2 space group. 
These symmetries arise in a variety of physical settings. 
For example, both symmetries are natural in systems of 
bosons, including situations where electrons form suffi¬ 
ciently tightly bound Cooper pairs. The former symme¬ 
try can arise in a Heisenberg or XY spin system if one 
ignores time reversal symmetry. The latter symmetry oc¬ 
curs in a Heisenberg model in a Zeeman field; the field 
naively breaks time reversal, but preserves a combination 
of time reversal and a 7r spin rotation perpendicular to 
the field axis. We focus on situations where G constrains 
the symmetry fluxes to be bosons, which simplifies the 
analysis; we show this occurs whenever time reversal or 
reflection symmetry is present. 

We do not discuss symmetries of the form G = U(l) x 
Zj x G spa ce. This important class of symmetries occurs in 
time-reversal symmetric XY or Heisenberg spin models. 
Application of our anomaly test for these symmetries is 
subtle (see Sec. VIII), and requires a more intricate anal¬ 
ysis that will be presented in a separate paper [40] . 

Partially for simplicity, and partially for its physical 
relevance, we concentrate on two-dimensional Z 2 topo¬ 
logical order, which means that the fusion and braid¬ 
ing of the anyon quasiparticles is the same as the de- 
confined phase of Z 2 gauge theory with gapped matter, 
or, equivalently, Kitaev’s toric code model [41]. SET 
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phases with Z 2 topological order are synonymous with 
gapped Z 2 quantum spin liquids (QSLs) [41-49], which 
are of current interest in part due to evidence that such 
a phase occurs in the S = 1/2 Heisenberg antiferromag- 
net on the kagome lattice [50-52]. While the symmetries 
we consider here are more relevant for other systems, Z 2 
QSLs can also occur in those systems. Showing that a 
given symmetry fractionalization pattern is anomalous 
constrains the possibilities for d = 2 Z 2 QSLs. 

As mentioned above, each anomalous symmetry frac¬ 
tionalization pattern we find provides a surface theory 
for a d = 3 bosonic TCI. Unlike the case of SPT phases 
protected by internal symmetry, there is not an existing 
theory of d = 3 bosonic TCIs, so it is particularly useful 
to obtain examples of such phases. We are able to ob¬ 
tain many such examples, and to discuss some of their 
physical properties, via their anomalous surface theories. 
It is not our goal to provide complete classifications of 
bosonic TCIs. 

For some bosonic TCIs, we can go beyond surface 
SET phases, and construct a dual vortex field theory 
for an anomalous surface superfluid. These superflu¬ 
ids, like some of the surface theories for bosonic topo¬ 
logical insulators studied in Ref. 35, are distinguished 
by non-trivial symmetry fractionalization of their vortex 
excitations [53-55]. The dual vortex field theories thus 
obtained are convenient to work with, and can be used 
to explore surface phase diagrams and phase transitions, 
which may be an interesting direction for future work. 

While it is not the focus of this paper, our approach 
can be used to study internal symmetries when G = 
U(l) xi G", and is complementary to existing approaches 
in that case. In particular, for G = U(l) x Z^, where Z^ 
is time reversal, our approach shows that certain fraction¬ 
alization patterns are anomalous, a result also obtained 
in previous works [35, 37]. The flux-fusion approach con¬ 
firms that result, without making assumptions about the 
form of the edge theory of d = 2 SET phases [37], or 
relying on a complete analysis of all possible phases of a 
surface field theory [35]. 

Table I summarizes the main results. Underlying the 
detailed results of the table are three distinct types of 
anomalies: 

1. For G = U(l) x Z[f, (T m ) 2 = —1 is anomalous, 
where T m gives the action of time reversal on vi- 
sons. 

2. Whenever G contains a U(l) x Zf subgroup, where 
Zf is reflection symmetry, (P m ) 2 = — 1 is anoma¬ 
lous, where P m gives the action of the reflection on 
visons. 

3. Whenever G contains a U(l) x Z^ subgroup, and 
also contains some discrete unitary operation g that 
commutes with the U(l) x Z^ subgroup, then 

T m g m = _ g m T m ( 1 ) 


is anomalous, where T m and g m give the action of 
T and g , respectively, on visons. For example, g 
can be a lattice translation or reflection. 

The first two types of anomalies can be understood in 
terms of dimensional reduction to d = 1 SPT phases, but 
it appears the third type of anomaly cannot be under¬ 
stood in this manner (Sec. V). 

We now provide some additional details in order to 
present Table I, followed by an outline of the remainder 
of the paper. As noted, we focus on Z 2 topological order, 
which supports four types of quasiparticle excitations, la¬ 
beled by 1, e, to, e. Of these, 1 particles are topologically 
trivial and can be created by local operators, while the re¬ 
maining particle types are anyons that cannot be locally 
created. We describe the fusion and braiding properties 
in Sec. Ill A. Here, we simply note that Z 2 topological 
order is realized in the deconfined phase of Z 2 gauge the¬ 
ory with gapped, bosonic matter, in which case e is the 
bosonic Z 2 gauge charge, to is the bosonic Z 2 gauge flux, 
and e is the fermionic charge-flux bound state. We will 
also refer to m particles as visons. 

We assume through the paper that symmetry does not 
permute the anyon species. In this case, the action of 
symmetry on the anyons is determined by giving the frac¬ 
tionalization class of e and to [56]. For each of e,m, the 
fractionalization class is an element of H 2 (G, Z 2 ). Here, 
this is specified uniquely by two pieces of information: 
1) whether the particle carries integer or half-odd integer 
U(l) charge, and 2) an element [w e ],[w m ] € H 2 (G', Z 2 ) 
that describes the action of G'. Each of e, to transforms 
as a projective representation of G' , and [w e ], [w m ] encode 
information about these projective representations that 
is a universal property of a SET phase (or surface SET 
phase). We always choose e to carry half-odd-integer 
charge, and to to carry integer charge. One motivation 
for this choice is that it describes most Z 2 QSLs that 
have been proposed to occur in fairly realistic models of 
spins or bosons. It can be shown, via a coupled layer 
construction [37], that all such symmetry fractionaliza¬ 
tion patterns (in fact, any consistent symmetry fraction¬ 
alization pattern) can occur as a surface of some d = 3 
SPT phase, which may be the trivial SPT phase (see 
Appendix F). Deciding whether the bulk SPT phase is 
non-trivial is equivalent to determining whether the cor¬ 
responding symmetry fractionalization pattern is anoma¬ 
lous. 

Under these assumptions, the flux-fusion anomaly test 
shows that some choices of [w m ] imply the symmetry frac¬ 
tionalization pattern is anomalous. This result is inde¬ 
pendent of [ uo e ], which does not play a role in the anomaly 
test. Column 2 of Table I is simply H 2 (G', Z 2 ), the set 
of all possible vison fractionalization classes for G' sym¬ 
metry. The anomaly test gives a subset of vison frac¬ 
tionalization classes that “test negative” for an anomaly 
and thus may occur strictly in d = 2. We refer to such 
classes as anomaly-negative ; they form a subgroup A f of 
H 2 (G', Z 2 ) given in column 3 of Table I. It is important 
to note that anomaly-negative fractionalization classes 


4 


may still be anomalous; the flux-fusion anomaly test can¬ 
not establish that a symmetry fractionalization pattern 
is non-anomalous. 

Finally, for a fixed [w e ], the anomaly test gives a set 
of distinct d = 3 SPT phases (one of which is always the 
trivial SPT phase), which are labeled by elements of the 
quotient S = H 2 (G', Z 2 )/Af, given in column 4 of the 
table. It is important to note that the anomaly test does 
not distinguish all SPT phases with a given symmetry, 
so column 4 does not give the full classification of such 
phases. 

We now give an outline of the remainder of the paper. 
Section II gives a simple, somewhat heuristic illustration 
of the anomaly test in the case of G = U(l) x (time 
reversal) symmetry. The anomaly test is then described 
in more detail and greater generality in Sec. III. First, 
Sec. Ill A describes the fusion and braiding properties 
both before and after gauging Z n C U(l). In Sec. IIIB, 
we describe the action of G' symmetry on the Z n flux fI, 
and use this to present the anomaly test. Especially for 
spatial symmetry, it is important for our analysis that f l 
is a boson, which is shown to be the case in Appendix B 
whenever time reversal or reflection symmetry is present. 

In Sec. IV, we apply the anomaly test to the examples 
of G = U(l) x Z l,G = U(l) x ZG = U(l) x pm, and 
G = (U(l) x Zj) x pi, where pm is a d = 2 space group 
containing translation and reflection operations, and pi 
is the d = 2 space group consisting only of translations. 
We find anomalous symmetry fractionalization patterns 
in each case. The first three of these symmetries have 
anomalies of type 1 and 2 as described above, which can 
be understood from the viewpoint of dimensional reduc¬ 
tion to d = 1 SPT phases, which is a different way to 
apply the anomaly test (Sec. V). In contrast, the last 
symmetry has type 3 anomalies that apparently cannot 
be understood in terms of dimensional reduction, as dis¬ 
cussed in Sec. V. 

Section VI describes how the results from the flux- 
fusion anomaly test can be used to identify and dis¬ 
tinguish some non-trivial d = 3 SPT phases, including 
d = 3 bosonic TCIs. As discussed in Sec. VII, some 
of the bosonic TCIs that we find can have an anomalous 
surface superfluid that preserves the G' symmetry. These 
anomalous superfluids are characterized by vortex excita¬ 
tions that transform projectively under the G' symmetry 
in a way that is not allowed strictly in d = 2. We describe 
how to construct dual vortex field theories that provide 
a convenient means to study the physical properties of 
these surface superfluids and neighboring surface phases. 

The paper concludes in Sec. VIII with a discussion 
of open issues raised by the present results. Some of 
the more technical aspects of our results are presented in 
several appendices, and, in Appendix G, the anomaly test 
is applied to a few more examples of symmetry groups. 


II. SIMPLE ILLUSTRATION OF THE 
ANOMALY TEST 

We begin by giving a somewhat heuristic illustration of 
the flux-fusion anomaly test, for the case of G = U(l) x 
Z^ 1 symmetry. This symmetry is chosen for simplicity, 
and for the fact that it has been previously studied using 
a different approach [37]. Here, we focus on conveying 
the intuition and some of the key ideas of our approach. 
A more rigorous and more general discussion follows in 
Sec. III. 

Here and throughout the paper, we assume d = 2 Z 2 
topological order, and that symmetry does not permute 
the anyon species. To specify the symmetry fractional¬ 
ization pattern, we need to give the fractionalization class 
for both e and m particles. For the present symmetry, we 
need to specify whether each particle carries integer or 
half-odd-integer U(l) charge, and whether it transforms 
as a Kramers singlet [(T a ) 2 = 1] or a Kramers doublet 
[(T a ) 2 = —1], where a = e,m, and T a gives the action 
of time reversal on anyon a. We denote particles with 
half-odd integer charge by C, and Kramers doublets by 
T, while 0 is used to indicate particles carrying trivial 
quantum numbers (integer charge and Kramers singlet). 
A fractionalization pattern is thus specified, for exam¬ 
ple, by the notation eCmT [37]; in this case, e particles 
carry half-odd-integer charge and are Kramers singlets, 
while m particles carry integer charge and are Kramers 
doublets. 

We restrict our attention to the case where e car¬ 
ries half-odd-integer charge and m carries integer charge, 
which includes four fractionalization patterns: eCmO, 
eCTmO , eCmT, and eCTmT. It is known that the for¬ 
mer two patterns are non-anomalous (can be realized in 
d = 2); this can be established, for example, via explicit 
construction of parton gauge theories. The latter two 
patterns were argued in Ref. 37 to be anomalous, via an 
approach that we contrast with ours at the end of this 
section. 

Our anomaly test is based on introducing fluxes fof 
the U(l) symmetry, where </> £ [0, 27 t). For the purposes 
of the present discussion, these fluxes are static point de¬ 
fects in space, obtained by modifying the Hamiltonian. 
The symmetry flux f2^ is defined by the following prop¬ 
erty: if Q is a local (i.e., non-anyon) excitation carrying 
unit U(l) charge, bringing Q counterclockwise around 
results in the statistical phase </>. We make the restriction 
0 < </> < 2-7T because and have the same mutual 

statistics with Q and thus carry the same symmetry flux. 

Given a fractionalization pattern, the flux-fusion 
anomaly test proceeds via two steps, which we summa¬ 
rize before proceeding. First, we study the fusion of sym¬ 
metry fluxes, and show that, roughly speaking, cf> = 2 tt 
flux is not trivial, but instead is a to particle excitation. 
Second, we consider the action of Z]f symmetry on sym¬ 
metry fluxes (f>, and ask whether it is possible to choose 
this symmetry action to be consistent with the assumed 
symmetry fractionalization of to, given the fusion prop- 
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erties of the fluxes. We will see there is an inconsistency 
if to is a Kramers doublet, so that eCmT and eCTmT 
are anomalous fractionalization patterns. 

First, to study the fusion properties of symmetry 
fluxes, we consider the mutual statistics of a flux f2^ 
with anyons e, e, to. We choose particular anyons e and 
e carrying U(l) charge 1/2, and to which is neutral un¬ 
der U(l). We could consider anyons with other allowed 
values of the charge (for example, there will also be e 
particles with charge —1/2), but this does not affect the 
results. Let be the statistical phase angle when 

anyon a is brought counterclockwise around the flux q i>. 
Then, given the assumed charge values for the anyons, 
we have 

©e.su = ©€,n 0 = | (2) 

©m,!l^ = 0. (3) 

To obtain some intuition for the fusion properties of 
the symmetry fluxes, suppose for the moment that we 
relax the restriction </> < 27 t. Then, if </> = 27r, we have 
formally 0 e ,n 2 „ = ©e,fi 27r = tt and & m ,n 2 „ = 0. Since 
I2 2w carries trivial symmetry flux (it has trivial mutual 
statistics with Q), it must be identified with one of the 
anyon quasiparticles. Putting <f> —> 2n in Eqs. (2) and (3), 
we have the identification f2 27r = to. Along the same 
lines, we can identify f = 1. 

We prefer to keep the restriction 0 < <j) < 2n, in which 
case essentially the same result can be obtained as fol¬ 
lows: Suppose that we have two n fluxes 12^. The total 
flux is 2-7T, which is equivalent to no symmetry flux at all. 
Therefore, we have the fusion rule 

= a , (4) 

where a is a quasiparticle excitation that carries no sym¬ 
metry flux, but may be a non-trivial anyon. The particle 
a can be identified by its mutual statistics with e, to and 
e, which follows from the additivity properties of statis¬ 
tics. For example, 

©e,a = ©e.SUSU — 2© e Q 7r = 7T. (5) 

Similarly, 0 ea = 7r and © mja = 0, which implies a = m 
and 

fln^-Tr = TO. (6) 

It should be noted that this result has a discrete char¬ 
acter and does not make use of the fact that U(l) is a 

continuous group. Indeed, the same result holds if we 
replace U(l) by the discrete group Z 2 . 

Next, we consider the action of time reversal symmetry 
T on the symmetry fluxes First, we observe that 
T does not change the value of the flux </>, because T 
commutes with U(l) rotations. Therefore, Q, v transforms 
either as a Kramers singlet or a Kramers doublet under 
time reversal. If we assume that to is a Kramers doublet, 
we now have a contradiction with Eq. (6): whether is 


a Kramers singlet or doublet, the composite must 

be a Kramers singlet. 

We have thus found that that eCmT and eCTmT are 
anomalous fractionalization patterns. This is true be¬ 
cause, in strict d = 2, it should always be possible to in¬ 
troduce U(l) symmetry fluxes and to view these as point 
objects, so the contradiction we obtained means that a 
fractionalization pattern cannot occur strictly in d = 2. 
On the other hand, on the surface of a d = 3 SPT phase, 
symmetry fluxes are line objects that penetrate into the 
bulk, and it may not be sensible to view them as point 
objects where they pierce the surface. Therefore, eCmT 
and eCTmT may occur on the surface of a d = 3 SPT 
phase. Indeed, this is the case, and was demonstrated in 
Ref. 37 via an elegant coupled layer construction. 

The above analysis is complementary to the approach 
of Ref. 37. There, among other results, Chern-Simons 
theory was used to construct chiral boson edge theo¬ 
ries for SET phases with Z 2 topological order and G = 
U(l) x Z 2 symmetry. For some symmetry fractional¬ 
ization patterns, including eCmT and eCTmT , it was 
shown that no corresponding edge theory can be con¬ 
structed, and it was concluded that these symmetry frac¬ 
tionalization patterns are anomalous. Strictly speaking, 
to draw this conclusion, one has to assume that the class 
of edge theories considered is in some sense sufficiently 
general, and, while this assumption seems reasonable, we 
do not know of an argument that this is the case. The 
flux-fusion approach requires no such assumption, and in 
the present case, its results agree with those of Ref. 37, for 
those fractionalization patterns where both approaches 
can be applied. 


III. FLUX-FUSION ANOMALY TEST: 

GENERAL DISCUSSION 

A. Gauging C U(l) symmetry 

The simple discussion of the anomaly test in Sec. II 
is based on inserting U(l) symmetry fluxes, which are 
static point defects in space. Because our objective is to 
consider crystalline symmetry, this approach is not ideal, 
because inserting a non-dynamical flux at some point in 
space will usually partially or fully break the crystalline 
symmetry. In addition, there is not an existing theory 
describing the action of G' symmetry on fluxes of the 
continuous U(l) symmetry. 

Therefore, we prefer to proceed by gauging a Z„ sub¬ 
group of the U(l) symmetry, for all integers n > 2. That 
is, we imagine minimally coupling our system to a dy¬ 
namical Z ra gauge held, where the Z n gauge group is iden¬ 
tified with Z n C U(l) global symmetry. In Appendix A, 
we give an explicit procedure showing that, for the sym¬ 
metry groups considered in this paper, it is possible to 
gauge this Z„ subgroup while preserving G'cG symme¬ 
try. The resulting theory is a gauged SET phase, where 
the symmetry flux behaves as a gapped, dynamical quasi- 
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particle excitation. This allows us to study symmetry 
fluxes without breaking crystalline symmetry. In addi¬ 
tion, we can build on existing results to describe the ac¬ 
tion of G' on the excitations of the gauged SET phase. 

We consider a d = 2 SET phase with Z 2 topological 
order and G = U(l) x G' symmetry. We now describe 
the fusion and braiding properties of the anyons of the 
SET phase. Fusion of anyons is described by the Abelian 
group A = Z 2 xZ 2l generated by e and to, which obey 
the relations 

e 2 = m 2 = 1 (7) 

e = era = me. (8) 

We assume that e carries half-odd-integer charge under 
U(l). Under Z n C U(l) symmetry, this means that 

{U e n T = - 1 , ( 9 ) 

where [7® is a unitary operator representing the action 
of a generator of Z n on a single e particle. Half-odd 
integer charge is only non-trivial for even n; if n is odd, 
then Eq. (9) can be trivialized by the allowed redefinition 
E/® —X —17®. Therefore, we restrict attention to even 
values of n. We also assume that m carries integer U(l) 
charge, so that under Z n we have (E/™)" = 1. The action 
of G' on e and m is characterized below in Sec. IIIB. 

To specify the statistics, we introduce some notation 
that will be particularly helpful in describing the gauged 
SET phase. For anyons a,b € A, let 9 a give the self¬ 
statistics angle of a, and let Q a ,b be the mutual statistics 
angle, where a is taken counterclockwise around b. These 
quantities satisfy the following general properties for any 
a, 6, c € A: 


o 

II 

e 

CD 

II 

(10) 

0a,a = 2 9 a 

(ii) 

O a , b = @6,a 

(12) 

®ab,c = ®a,c H - ®b,c 

(13) 

Oab = 0 a ~\~ 6b ~\~ O a,6- 

(14) 

These and other equations for 6 a and Q a ,b 

are always 

understood to be true modulo 27 t. The statistics of Z 2 
topological order is then fully specified by 

9 e =9 m = 0 

(15) 

®e,ra = 7T • 

(16) 


These equations say that e and to are bosons with 0 em = 
7r mutual statistics. 

We now consider the gauged SET phase, obtained by 
gauging Z n C U(l). The anyons of the gauged SET phase 
are Abelian; this follows from Eq. 399 and the surround¬ 
ing discussion of Ref. 32. The fusion rules are described 
by the Abelian group Ag, which is generated by e, to, Q, 
and f 1. Here, Q is the unit Z„ symmetry charge, which 
is a local excitation of the un-gauged theory, but is now 
an anyon in the gauged SET phase. Q is the unit Z„ 
symmetry flux. Upon gauging Z n , the e and to sectors in 


the un-gauged theory each break into n different sectors 
with distinct Z n symmetry charge. In the gauged SET 
phase, e and to each correspond to a particular choice 
among such subsectors. The choice of subsector is arbi¬ 
trary, and can be changed by redefining e or to by binding 
symmetry charges; for example e —> Qe is an allowed re¬ 
definition. There is also arbitrariness in the choice of 
symmetry flux, which can be redefined by Q —> Qfl, or 
by Q —> afl, where a is an anyon of the un-gauged theory. 

The fusion rules are 


Q n = 1 

(17) 

e 2 = Q 

(18) 

m 2 = 1 

(19) 

= aQ k . 

(20) 


Equation (17) is obvious. Equations (18) and (19) corre¬ 
spond to making a particular choice of e and to among 
the possible subsectors. The most important fusion rule 
in our analysis is Eq. (20). There, a is an anyon of the 
un-gauged theory to be determined, and k is some as yet 
unknown integer satisfying 0 < k < n. This equation 
expresses the fact that Q n carries no Z n symmetry flux, 
but otherwise, at this stage in the analysis, could be an 
arbitrary particle in the gauged SET phase. 

In order to fix the fusion rule Eq. (20), we consider the 
statistics of the gauged SET phase. We have 


0 e =0 m = 0 ( 21 ) 

0e,m = 7T (22) 

0 Q = &e,Q = © m ,Q = 0 (23) 

27T 

0Q,o = - (24) 

n 

0e,n = - + PeTT (25) 

n 

0m,f2 = p rn 7T. (26) 


Here, Eqs. (21) and (22) are the braiding statistics for the 
un-gauged SET phase. Equation (23) holds because the 
symmetry charge Q must have trivial braiding with itself 
and with anyons of the un-gauged theory. Equation (24) 
is the defining property of the symmetry flux fl. Finally, 
Eqs. (25) and (26) follow from Eqs. (18) and (19), re¬ 
spectively, with unknown parameters p e ,Pm = 0,1. 

We redefine e and to to set p e = p m = 0. For example, 
if p e = 1, we redefine e —X Q”/ 2 e. This leaves the fusion 
rules unchanged, and results in 0 Ej n = 7r/n, without 
modifying the other statistics angles. 

Now, we use the statistics to constrain the flux fusion 
rule, Eq. (20). Using 0 mj n = 0, we have 0 m ,n™ = 0. 
Consistency with Eq. (20) then requires either a = 1 or 
a = to. Similarly, Q e ,n = 7r jn implies 0 ei n™ = vr, which 
requires either a = to or a = e. Therefore, a = to, and 

n n =mQ k . (27) 

So far, we have not mentioned 9q, the self-statistics of 
the symmetry flux. Unlike the other statistics angles, this 
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parameter does not follow immediately from our assump¬ 
tions, but it can be related to the integer k appearing in 
Eq. (27). First, Eq. (27) implies that Q n is a boson, so 
0q™ = n 2 8 n = 0 , and therefore 


On = 



(28) 


for some integer q satisfying 0 < q < n 2 . In fact, we can 
further restrict the range of q. To see this, we make the 
redefinition Q — > efl and m —> Q n ^ 2 m, which preserves 
the fusion rules, and leaves all the statistics angles un¬ 
changed except On- The effect of this redefinition is to 
shift q —>• q+n/ 2 , which allows us to restrict 0 < q < n/2. 

We can now relate q and k by noting that = 

2 uOq = Anq/n, and also ©n,£ 2 " = ©n,mQ fc = 27r k/n, so 
that Airq/n = rk/n. This has no solution if k is odd, so 
k must be even. Given the restrictions on the range of k 
and q, the unique solution for q is then q = k/ 2 , and we 
have shown 


and the statistics are specified by 


0e = e m = 0 
©e,m = 7T 

0 q = 0 fi .o = 8m.O = 0 


y e,Q 


m,Q 


<=> 2lT 
°Q,Q - - 

Be.O = - 

n 

©m,n = 0 
0 o = O. 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 


These are precisely the fusion rules and statistics of Z 2 „ 
gauge theory, or, equivalently, the Z 2 ra version of the toric 
code model. For Abelian anyons, fusion rules and statis¬ 
tics are enough to uniquely specify the unitary modular 
tensor category that describes a theory of anyons [57, 58] . 
Therefore, the theory of anyons in the gauged SET phase 
is identical to that in the Z 2n toric code. 


B. Symmetry action on m, Q and the anomaly test 


where k is even and satisfies 0 < k < n. In particular, 
for n = 2 we have Q 2 = m, as stated in Sec. II. 

Physically, we expect k to parametrize the quantized 
Hall response. Inserting 27 t flux at some point in space 
produces a local charge accumulation of a xyi in appro¬ 
priate units. If we view fusion of n fluxes Q as equivalent 
to a dynamical process where n fluxes are inserted, then, 
because m is neutral under Z n , Eq. (27) implies 

k = <J xy mod n. (30) 

This physical interpretation of k leads us to expect k = 0 
whenever G' symmetry forbids a quantized Hall response. 
Indeed, in Appendix B, we show that k = 0 whenever G' 
contains time reversal or spatial reflection symmetry. 

Whenever k = 0, by Eq. (29), fl is a boson. This will 
enable a simple description of the action of G' symme¬ 
try on f! and in , so from now on we will always assume 
conditions are such that we can take fl to be a boson. Un¬ 
der this assumption, we collect here the properties of the 
gauged SET phase obtained from the discussion above. 
The fusion rules are 


Q n = 1 

(31) 

e 2 = Q 

(32) 

TO 2 = 1 

(33) 

n n = m, 

(34) 


In order to apply the anomaly test, we first have to 
characterize the action of G' symmetry on the anyons 
of the un-gauged SET phase [19]. In general, the frac- 
tionalization class of e or m is an element of the group 
H 2 (G, Z 2 ). In the present case, as is shown in Ap¬ 
pendix C, it is enough to specify separately the action 
of U(l) and G' on each of e and m. That is, there is no 
additional information associated with interplay between 
U(l) and G'. 

Each of e, m transforms under a Z 2 projective repre¬ 
sentation of G' denoted F e , F m , respectively. We focus 
on m particles; the corresponding equations hold for e 
particles with trivial modifications. For 51,52 £ G 1 , we 
have 

r m (5i)r m (52) = cc m ( 5 1 , 5 2 )r m ( 5 1 5 2 ), (42) 

where w m (gi,g 2 ) £ Z 2 is called a Z 2 factor set. The 
corresponding object for e particles is denoted u> e . Asso¬ 
ciative multiplication of the r m ’s implies 

w m (5i,5 2 )w m (5i52,53) = w m (5i,525 3 )w m (52,53)- (43) 

In general, any Z 2 -valued function 1 ^( 51 , 52 ) satisfying 
Eq. (43) is called a Z 2 factor set. 

Physical properties are unchanged under a redefinition 
r m ( 5 ) —► A” 1 ( 5 )r m ( 5 ) for A( 5 ) £ Z 2 , which induces a 
projective transformation on the factor set, 

w m (5i,5 2 ) -> A - 1 (5i)A - 1 (5 2 )A(5i5 2 )w m (5i,5 2 ). (44) 

Here, A - 1 ( 5 ) = A( 5 ), but the inverse signs are kept to 
expose the formal similarities with the discussion of sym¬ 
metry action on H, below. Equivalence classes of factor 
sets under such projective transformations are denoted 
[w m ]z 2 , and are the distinct fractionalization classes of 


to. The Z 2 subscript reminds us that both w m and 
the projective transformations A take values in Z 2 . In 
the language of group cohomology theory, fractionaliza- 
tion classes [w m ]z 2 are elements of the Abelian group 
H 2 {G', Z 2 ), the second group cohomology of G' with 
Z 2 coefficients. The group multiplication in iJ 2 (G',Z 2 ) 
is obtained from multiplication of functions; that is, if 
Wab(gi,g 2 ) = Va(gi,g 2 )ub(gi,g 2 ), then [w a ]z 2 [w h ]z 2 = 
[w a b] Z 2 - 

Considering all symmetries together, the symmetry 
fractionalization pattern of the SET phase can be de¬ 
noted eC[w e ]TO.0[w m ], where C (0) indicates that e (to) 
carries half-odd-integer (integer) U(l) charge. When us¬ 
ing this notation, to avoid cumbersome expressions, we 
drop the Z 2 subscript for the fractionalization classes. 

The flux-fusion anomaly test will be able to determine 
that eC[ca e ]mO[u} m ] is anomalous for certain choices of 
[w m ]z 2 , independent of [w e ]z 2 - When the anomaly test 
does not find an anomaly, we say that a symmetry frac¬ 
tionalization pattern is anomaly-negative. This termi¬ 
nology recognizes that the flux-fusion anomaly test is 
not expected to detect all possible anomalies, and some 
anomaly-negative fractionalization patterns can still be 
anomalous. 

To proceed, we now consider the gauged SET phase, 
and characterize the action of G' symmetry on Cl. First, 
we need to consider the possibility that some operations 
may permute the anyons of the gauged SET phase, and, 
in particular, may map Cl to some other anyon. For some 
operation g G G', let g*Cl denote the anyon resulting from 
applying g to Cl. If g commutes with U(l), is unitary, and 
is either an internal symmetry or a proper space group 
operation, then g-kCl = Cl. This follows from the fact that 
such an operation leaves Q , e and to invariant, and also 
leaves the statistics invariant; that is, 0 g * a , g *b = O a ,b- 
However, it is not the case that all g G G' leave Cl in¬ 
variant; in particular, we will be interested in time rever¬ 
sal and reflection symmetry. These operations may send 
Cl or Cl 1 —> fl 2n_1 , depending on whether the opera¬ 
tion in question commutes with U(l), as is discussed in 
detail in Appendix B. 

Because some operations in G' may not preserve the 
anyon type of Cl, in describing the action of symme¬ 
try, we have to go somewhat beyond the framework de¬ 
veloped in Ref. 19. We introduce field operators 
(jfe = 1,..., 2 n — 1). Each xjjf. is a many-component ob¬ 
ject, with components not explicitly written, where each 
component creates a Cl k particle in some state. In partic¬ 
ular, il> n creates a to particle. These field operators are 
non-local objects. However, because all the Cl k particles 
are bosons and have bosonic mutual statistics, the non¬ 
local character of i/ife is not expected to play a role in the 
following discussion. It is also convenient to collect all 
the field operators into the object 4/ = (tpi • • • ip 2 n- i). 

All physical states and local operators are invariant un¬ 
der Z 2 n gauge transformations implemented by the uni¬ 
tary operator C?[A], for A G Z 2ra , which acts on the field 


operators by 


G[M*PkG[ A]” 1 = A Vfc- (45) 

For a symmetry operation g G G', we denote the corre¬ 
sponding unitary or anti-unitary operator by S(g ), which 
acts on field operators by 

g : * Sig^Sig)- 1 . (46) 

The operators S(g) form a representation of G' up to Z 2ra 
gauge transformations, that is 

S(gi)S(g 2 ) =Q[4>n{gi,g2)]S(gig2), (47) 

for ^ 1 ( 31 , 32 ) G Z 2 n- This is the most general multipli¬ 
cation law consistent with the requirement that S(g) act 
linearly on local operators, for example (schematically), 
(V’l) 2 "- Mathematically, we have defined a kind of gen¬ 
eralized projective representation, which is similar to but 
not identical to the projective representation describing 
the action of symmetry on m [Eq. (42)]. 

The crucial difference between S(g) and more familiar 
projective representations is that, in general, S(g) does 
not commute with the gauge transformation C?[A]. We 
note that some symmetries g G G' map g : Cl >—>• d 2n ~ 1 . 
We keep track of this information by defining 


s(g) 


“hi, q \ 1 —y 

— 1 , q : Q 1 —y Sl 2n_1 


(48) 


In addition, some operations in G' may be anti-unitary, 
so we define 


u(g) = 1; 5 Unitary . . (49) 

1—1, g anti-unitary 

We note that both s and u are group homomorphisms 
mapping G' —> Z 2 . We then introduce the function 

t(g) = s(g}u{g). (50) 

By considering the action of S ( g ) and gauge transforma¬ 
tions on field operators, it is straightforward to show 

S(ff)G[A] = S[A t(9) ]S( 5 ), (51) 

which shows that t{g) characterizes the non¬ 

commutativity of S(g) and gauge transformations. 
We thus refer to S(g) as a f-twisted Z 2n projective 
representation of G'. 

Equation (51) allows us to use associativity of the 
product S(gi)S(g 2 )S(g 3 ) to derive the associativity con¬ 
dition on <j> n , 

< t > n { gi , 52)^(5152, 53 ) = <Msi, 3233)^(32, 33 )F (si) - 

(52) 

We refer to <fi n , and, indeed, any Z 2 „-valued function sat¬ 
isfying Eq. (52), as a t- twisted Z 2ra factor set. Parallel¬ 
ing the discussion of ordinary projective representations 
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above, we are free to redefine S(g) by a gauge transfor¬ 
mation, 

S(g) -»• GlX-'mSig). ( 53 ) 

This induces a projective transformation on the factor 
set, 

(t>n(gi,g2) -> A _1 (5i)[A(5f2)] _t(9l) [A(5i32)]^n(5i,52)- 

(54) 

Equivalence classes [<f>n]z 2ll of factor sets under such 
transformations characterize the action of G' symmetry 
on f2. These equivalence classes are elements of the coho¬ 
mology group Z 2n ), where the t subscript denotes 

the non-trivial action of G' on the Z 2 n coefficients, en¬ 
coded in the function t(g). We refer to this as f-twisted 
cohomology. We note that, for G' finite, on-site and uni¬ 
tary, we have recovered a special case of the twisted coho¬ 
mology theory used to describe the action of symmetry 
on anyons in the category-theoretic description of SET 
phases [32, 59]. 

In fact, [4> n \ z 2 „ simultaneously characterizes the action 
of G' on all particles f2 fc obtained by fusing fl’s together. 
This includes m = f l n . The action of symmetry on to is 
given by considering the action of S(g) on ip n -, in partic¬ 
ular, 

S(gi ) S {g 2 ) ijj n s (c/2 ) _11 S(gi ) “ 1 

= [<t>n{gi, g 2 )] n S(gig 2 )ipnS{gig 2 )~ 1 
= Vm(gi, g2)S(gig 2 )4>nS{gig 2 )~ 1 ■ (55) 

Therefore, we have shown 

(51,02) = [0n(0i,02)] n - (56) 

Equation (56), which holds for all even n > 2, is the 
crucial equation underlying the anomaly test. The essen¬ 
tial idea is to take advantage of the fact that 12 is a “nth 
root” of to in the gauged SET phase, and to ask whether 
a given symmetry action on to can be consistently ex¬ 
tended to a symmetry action on its nth root 12. If not, 
then an anomaly has been detected. 

In more detail, the logic is as follows: Given [w m ]z 2 , 
we choose some particular factor set u! m (gi, g 2 ) in the de¬ 
sired equivalence class (the particular choice within the 
class does not matter). Then, for each even n > 2, we 
ask whether it is possible to solve Eq. (56) for 4 > n {g\,g 2 ), 
where <p n is required to satisfy Eq. (52). If for any even 
n > 2 , a solution fails to exist, the symmetry fraction- 
alization pattern is anomalous. If a solution exists for 
all even n > 2 , the symmetry fractionalization pattern is 
anomaly-negative. 

Equation (56) immediately implies that anomaly¬ 
negative to particle fractionalization classes form a sub¬ 
group that we denote TV C H 2 (G Z 2 ). 

At first glance, it might appear cumbersome to apply 
the flux-fusion anomaly test. Fortunately, it is not nec¬ 
essary to consider Eq. (56) directly for every even n > 2. 
Instead, there is a simple and easily computable char¬ 
acterization of which [w m ] Z2 are anomaly-negative. To 


describe this characterization, we first note that ui m can 
be viewed as a f-twisted U(l) factor set. This means that, 
given w m (< 7 i, g 2 ), we allow for projective transformations 

w m (si,g 2 ) A ~ 1 (gi)[\{g 2 )]- t< ' 9l) [\{gig 2 )\u m (gi,g 2 ), 

( 57 ) 

where A(g) £ U(l). The corresponding equivalence class 
under these transformations is denoted [w m ]u(i), and is 
an element of the cohomology group H 2 (G\ U(l)). For¬ 
mally, there is a map p 2 : H 2 (G\ Z 2 ) —> H 2 (G\ U(l)) 
defined by p 2 ([w m ]z 2 ) = [w m ]u(i)- (In Appendix D, it is 
shown that p 2 is well-defined, is a group homomorphism, 
and is unique in a certain natural sense.) 

Intuitively, it seems natural for cohomology with U(l) 
coefficients to arise out of the flux-fusion anomaly test. 
Ultimately, it ought to be possible to dispense with gaug¬ 
ing Z n C U(l) for all even n , in favor of working directly 
with continuous U(l) fluxes. Either approach should give 
the same results, so we speculate that the H 2 (G\ U(l)) 
cohomology may describe the action of G’ symmetry on 
U(l) fluxes. Moreover, as discussed in more detail in 
Sec. VII, [w m ]jj(ij does have a nice physical interpreta¬ 
tion: it characterizes the symmetry fractionalization of 
vortex excitations in a superfluid. This allows us to ob¬ 
tain results on anomalous d = 2 superfluids. We note 
that f-twisted U(l) cohomology also appears in the coho¬ 
mology approach to SPT phases with time reversal sym¬ 
metry, where anti-unitary operations act non-trivially on 
the U(l) coefficients, and the cohomology groups are de¬ 
noted by H n (G, U T (1)) [16]. 

Anomaly-negative to particle fractionalization classes 
[w m ]z 2 are fully characterized by the following theorem, 
which is proved in Appendix D. 

Theorem 1. If H 2 (G\ U(l)) = U(l) fc x A, where A is 
a finite product of finite cyclic factors, then the symme¬ 
try fractionalization pattern eC[w e ]TO 0 [a; m ] is anomaly- 
negative if and only if [w m ]u(i) = p 2 ([w m ]z 2 ) lies in the 
connected component of Hf(G', U(l)) that contains the 
identity element. 

The assumption on the form of Hf{G', U(l)) is true 
for all the examples we have considered, and we believe 
it is likely to be true in general. 

This theorem allows us to apply the flux-fusion 
anomaly test via the following procedure: 

1. Compute the group H 2 (G', Z 2 ) of to particle frac¬ 
tionalization classes under G' symmetry. Find a 
convenient explicit parametrization of distinct frac¬ 
tionalization classes [w m ]z 2 - 

2. Compute the f-twisted cohomology group 
Hf (G', U(l)), and find an explicit parametrization. 

3. Find the map p 2 discussed above, for which 
[w m ]u(l) = P 2 ([^m]z 2 ) ■ 

4. For each to particle fractionalization class [w m ]z 2 , 
determine whether [cc TO ]u(i) can be continuously 
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deformed to the identity element of i7 2 (G", U(l)). 
If this is impossible, the fractionalization pattern 
eC[w e ]mO[u; m ] is anomalous, for any [w e ]z 2 . 

5. The results for a given symmetry G = U(l) x G' 
can be summarized by describing the m particle 
fractionalization classes for which eC[w e ]mO[w m ] is 
anomalous. 

This procedure is illustrated in detail, and made more 
concrete, in the examples presented in Section IV and 
Appendix G. 

IV. EXAMPLES 

In this Section, we apply the flux-fusion anomaly test 
in a few cases, in order of increasing complexity. In each 
case, we fix a symmetry G , and follow the procedure 
outlined in Sec. IIIB. A crucial aspect is the calculation 
of second cohomology groups, for G' presented in terms 
of generators and relations. We illustrate our approach to 
these calculations in each example, leaving a more careful 
mathematical justification to Appendix E. 

These examples enable a more concrete discussion of 
bosonic TCIs in Sec. VI, and anomalous d = 2 superfluids 
in Sec. VII. A number of other examples are considered 
in Appendix G. 

A. G = U(1)xZT 

We begin with the case of G = U(l) x Z]f symme¬ 
try, that was already discussed in Sec. II and in previous 
works [35, 37]. This symmetry is simple enough to an¬ 
alyze using Eq. (56) directly; that approach, in fact, is 
essentially identical to the treatment in Sec. II. However, 
to pave the way for more complex examples, we follow 
the procedure outlined in Sec. IIIB. 

It is convenient to present the group G' = Z]T in terms 
of generators and relations. Here, this is trivial; the single 
generator T obeys the relation T 2 = 1. Next, we consider 
a general Z 2 projective representation giving the action 
of G' on a m particle. The generator is now written T m , 
and the relation becomes 

(T m ) 2 = ( 58 ) 

for cr™ £ Z 2 . We are allowed to redefine T m —> —T m , 
but this does not affect cr™. Therefore, because cr™ = ±1 
is invariant under projective transformations, we can ten¬ 
tatively conclude that it labels two distinct fractionaliza¬ 
tion classes [ui m ]z 2 . To be sure this conclusion is cor¬ 
rect, we need to check that each choice of a™ in fact 
corresponds to a factor set w m (< 7 i, g 2 ), for gi,g 2 £ Z^. 
This can be accomplished by exhibiting a projective 
representation for each choice of cr™. In the present 
case, these representations are just familiar Kramers sin¬ 
glets (cr™ = 1 ) and doublets (cr™ = — 1 ). Therefore, 


H 2 (G\ Z 2 ) = Z 2 , with cr™ = ±1 explicitly parametrizing 
the cohomology group, and labeling the distinct fraction¬ 
alization classes [w m ]z 2 ■ 

Next, we have to compute Hf(G r , U(l)). To do this, 
we consider a general t-twisted U(l) projective represen¬ 
tation of G' = Z%, again in terms of generators and rela¬ 
tions. We denote the generator by T f . We also have to 
specify the function f(g); it is sufficient to give the values 
of t for the generators, and in this case, t(T) = —1. The 
relation becomes 

(T*) 2 = a T . (59) 

Here, ax is short-hand for the gauge transformation 
G\oit\ , with ax £ U(l). So, for example, we can write 
T 4 Q!t = ot^ 1 T t ■ It is important to note that we can ad¬ 
just the phase of the generator, by redefining T* —> AT*, 
but this leaves ax unchanged. 

Because ax is invariant under projective transforma¬ 
tions, it is tempting to conclude that ax £ U(l) labels 
distinct equivalence classes [w]u(i) G -ff 2 (G', U(l)). How¬ 
ever, this conclusion is not correct, because the possible 
values of ax are constrained. That is, there does not ex¬ 
ist a t-twisted U(l) factor set for arbitrary ax £ U(l). 
To see this, we conjugate both sides of Eq. (59) by T*, 
and readily obtain (T *) 2 = a^ 1 . This is consistent only 
if ax £ Z 2 . 

As before, we need to verify that both choices ax = 
±1 actually give rise to t-twisted U(l) factor sets. The 
same Kramers singlet and doublet representations can 
be viewed as t-twisted U(l) projective representations, 
so, once again, we can exhibit a representation realizing 
each choice of ax- Therefore, H 2 (G', U(l)) = Z 2 , which 
is explicitly parametrized by ax £ Z 2 . 

To find the map p 2 giving [w m ]u(i) in terms of [w m ]z 2 , 
suppose we have a Z 2 projective representation as de¬ 
scribed in Eq. (58) with some value of a™. This Z 2 
projective representation can immediately be viewed as 
a t-twisted U(l) projective representation, with ax = 
P2(&x) = a T- I n ^is case, then, p 2 : Z 2 —>• Z 2 is the 
identity map; more non-trivial examples will arise for 
other symmetries. 

To conclude, we see that cr™ = 1 is anomaly¬ 
negative, because p 2 (l) = 1- On the other hand, 
cr™ = — 1 is anomalous, because p 2 (— 1 ) = — 1 , which 
is not continuously connected to the identity element in 
H 2 (G\ U(l)) = Z 2 . The group A f of anomaly-negative 
vison fractionalization classes is thus trivial, Af = 7L\. 

B. G = U(l)xZf 

Next, we consider the case of a single lattice reflec¬ 
tion symmetry (Zf), which commutes with the U(l). Of 
course, any realistic system with reflection symmetry will 
also have a larger space group, including translation sym¬ 
metry. A physically reasonable viewpoint is to imagine 
that we are interested in a system that has additional 
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space group symmetry beyond Zf, but we are “forget¬ 
ting” about the rest of the symmetry, and only making 
use of a U(l) x Z^ subgroup in our analysis. 

Our discussion parallels the treatment given above for 
time reversal symmetry. The group G' = Zf is generated 
by P, which obeys the relation P 2 = 1. Acting on m 
particles, the generator is written P m and obeys 

(P m ) 2 = (jp , (60) 

with (jp € Z 2 . As before, < 7 p is invariant under P m —» 
—P m . Both choices of <jp can be realized; for example, 
we can choose one-dimensional representations P m = 1 
(for <jp = 1) and P m = i (for crj? = -1). Therefore, 
H 2 {G', Z 2 ) = Z 2 , parametrized by cr™. A physical con¬ 
sequence of dp is that, if a pair of m particles are created 
and moved to reflection symmetric points, the resulting 
wave function has a reflection eigenvalue of <Tp, relative 
to the reflection eigenvalue of the ground state. 

Next, we consider a general t-twisted U(l) projective 
representation generated by P*, obeying the relation 

(P ‘) 2 = a P , (61) 

with ap £ U(l). In Appendix B, it is shown that P maps 
the symmetry flux tt to since P is unitary, this 

implies f(P) = —1. 

At this point, the analysis is mathematically identical 
to that given for G = U(l) x Z^ symmetry in Sec. IV A. 
That is, we have ap £ Z 2 , and ff 2 (G', U(l)) = Z 2 . In 
addition, ap = p 2 (erp) = a™. Therefore, a™ = 1 is 
anomaly-negative, while er™ = — 1 is anomalous, and the 
group of anomaly-negative vison fractionalization classes 
is Af = Zi. Introducing notation similar to that used in 
Sec. II for time reversal symmetry, we have found that the 
symmetry fractionalization patterns eCmP and eCPmP 
are anomalous, where P denotes an anyon for which P 2 = 
- 1 . 


C. G = U(l) x pm 

We now move on to an example in which G' is a d = 2 
space group. We choose G' = pm , which is generated by 
translations T x , T y , T" 1 , T” 1 , and a reflection P x \ these 
operations are illustrated in Fig. 1, and obey the relations 


T x T y T~X X = 1 

(62) 

TyP X T~ 1 P X = 1 

(63) 

P 2 = 1 

(64) 

T X P X T X P X = 1 , 

(65) 


which define the group pm. 

Acting on m particles, the general form of the relations 
is 


rpm rprn rpm—1 rpm — 1 
x -L y x y 

_ m 

®txty 

( 66 ) 

rpm pymrpm— 1 pm 

^ y r x y *-x 

_ m 

u typx 

(67) 

{PX 

_ m 

®px 

( 68 ) 

rpm p>m rpm pm 

_ m 

( - f txpx’> 

(69) 



FIG. 1. Illustration of the operations generating the d = 2 
space group pm.. A square lattice, which is invariant under 
pm symmetry, is shown to aid visualization. T x and T y are 
translations by one lattice constant along the x- and y- axes, 
respectively. P x is a reflection, with axis indicated by the 
left-hand vertical dashed line. The group pm has two types 
of reflection axes, with P x being of one type, and T X P X being 
of the other type. The reflection axis for T X P X is shown as 
the right-hand vertical dashed line. 

with the cr m ’s taking values in Z 2 . All the cr m ’s are invari¬ 
ant under projective transformations of the generators 
( e.g. T x —» — T x ), which suggests P 2 (G',Z 2 ) = 7L\. Be¬ 
cause these relations are a subset of those used to present 
the square lattice space group in Ref. 19, it follows from 
Appendix A of that work that all 2 4 possible choices of 
the cr m ’s indeed correspond to a factor set w m , and in¬ 
deed H 2 {G\ Z 2 ) = Z| (see also Appendix G 1 of this pa¬ 
per). The fractionalization classes are thus parametrized 

v.. r. , 1 _ / _m _m _m _m \ 

[uJm JZ2 \ u txtyi u typxi u pxi u txpx) m 

Now we need to compute H 2 {G\ U(l)), noting that 
t{P x ) = —1, while t{T x ) = t(T y ) = 1. The general form 
of the relations in a f-twisted projective representation is 


TXTtX- 1 = atxf y ( 7 °) 

TXX~ lp * = a typx ( 71 ) 

{Px ) 2 — a p X (72) 

Tt x P i T l pt x = a txpx, (73) 


where the a’s take values in U(l). If we redefine Ty —>■ 
{oitypxY^Ty, this sends at ypx —► 1 , leaving the other a’s 
unchanged. In addition, the other a’s are unchanged by 
redefinition of the other generators, so we have arrived at 
a convenient canonical gauge choice to describe a general 
t-twisted projective representation. 

Next, conjugating Eq. (72) by P X1 we find a px £ {±1}, 
and similarly find at xpx £ {±1} by conjugating Eq. (73) 
by T*P f x . This suggests that H 2 {G’, U(l)) = U(l) x 
Z 2 x Z 2 , with the elements of the cohomology group 
parametrized by [w] u(1 ) = {a txty 5 Qpxi &txpx) • 

To verify this, we need to exhibit t-twisted projective 
representations that correspond to a generating set of 
H 2 {G', U(l)). (It is enough to exhibit a generating set, 
because the corresponding factor sets can then be mul¬ 
tiplied to obtain a factor set with arbitrary cohomology 
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class.) We introduce field annihilation (creation) opera¬ 
tors v r (uj) for some fictitious particles residing on the 


sites r = (x,y) of the square lattice. Each v r is a two- 
component vector. The generators are chosen to act on 
the field operators by 

T X V r T x = { atxty ) v 9tx^r+x (74) 

TyVrTy = gtyV r -\-y ( 75 ) 

PxV r Px 1 = 9pxVp xr , (76) 

where a txty G U(l), P x r = (-x,y), and gt x ,9ty,g P x are 
2x2 unitary matrices. Gauge transformations act on the 
field operators by 

Q[\\v r g- x [\]=e iXrT ‘v r (77) 

g[\}vlg- 1 ix\ = e- iX ° z vl (78) 


for A G U(l), with a z one of the 2x2 Pauli ma¬ 
trices, and where u) denotes Hermitian conjugation of 
the components of v r , but does not include transpo¬ 
sition in the two-component space. Choosing g tx = 
9ty = g P x — 1 gives a continuous family of represen¬ 
tations with [w] u(1) = {a tx ty, 1,1)- Next, a txtv = 1, 
gt y = i, g px = icr v , gt x = <r 2 is a representation with 
Mu(i) = (1, —1,1). (Here, again, a x ,a v ,a z are the 
usual 2x2 Pauli matrices.) Finally, atxty = 1, gty = 1, 
g P x = cr x , gtx = c 2 has Hu ( i) = (1,1, —1)- The fac¬ 
tor sets of these three families of representations are a 
generating set for H?(G', U(l)) = U(l) x Z 2 x Z 2 . 

To find the map p 2 , we begin with [w m ]z 2 = 
{vZty^tZpx^Px’VtZpx)- Viewing a corresponding pro¬ 
jective representation as a t-twisted U(l) projective rep¬ 
resentation, we can redefine T™ —► (o"J^ ) 1 // 2 T™, which 

sets —> 1 , and thus puts this projective representa¬ 

tion in the canonical gauge described above. Therefore, 
we have found 

{atxty > a px , atxpx) = P2([ t a m ]z 2 ) = { <J txtyl <7 pxl Cy tx P x)- 

. (79) 

This can be continuously deformed to the identity in 
H?{G', U(l)) if and only if < 7 ™ = a% px = 1. 

Therefore, we have found that the fractional- 
ization pattern eC[w e ]rnO[w m ], with [w m ]z 2 = 

is anomalous if er™ = -1 
or cr^ px = —1 (or both). Equivalently, we can 

observe that anomaly-negative [w m ]z 2 are given by 
[w m ]z 2 = (v^ty^ujpx, 1 ,!), which form the subgroup 
Afc ff 2 (G',Z 2 ), with A7= (Z 2 ) 2 . 

D. G = (U(l) x 1%) x pi 

The group pi is the d = 2 space group consisting only 
of translation symmetry. Here, we consider this sym¬ 
metry combined with time-reversal, which enters via the 
semi-direct product U(l) x Z|\ This example is straight¬ 
forward to analyze by following the steps in Sec. IV C 
and Appendix G 1, so we only quote the results. 


We have G' = pi x Z^, which is generated by trans¬ 
lations T Xl T y , their inverses, and time reversal 77 The 
relations are 


TxTyT-'T- 1 = 1 

(80) 

r 2 = 1 

(81) 

TT X = T X T 

(82) 

TTy = TyT- 

(83) 

The m particle symmetry fractionalization is 

specified by 

rprnrprnrprn—lrprn— 1 _ m 

x y x y u txty 

(84) 

( / ) = cr T 

(85) 

m rpm^rm 
1 1 x ~ <J Ttx 1 x 1 

(86) 

'-pmrpm m rpm'-pm 

1 1 y ~ <J Tty 1 y ’ 

(87) 


where the er m ’s are Z 2 -valued phase factors. The group 
of vison fractionalization classes is H 2 {G', Z 2 ) = (Z 2 ) 4 . 
The anomaly-negative fractionalization classes are those 
with <J™ tx = <Jpty = V and form the group M = (Z 2 ) 2 . 
We thus have an anomalous fractionalization pattern 
when (Jp tx = —1, <r™ t = — 1, or both. The disjoint 
sets of SPT phases distinguished by the anomaly test are 
labeled by elements of S = H 2 {G', Z 2 )/A7 = (Z 2 ) 2 . 

It is interesting to note that, in this case, we find 
anomalies involving time reversal that cannot be under¬ 
stood in terms of U(l) x Z^ subgroups of G; these are 
type 3 anomalies as discussed in Sec. I. These anoma¬ 
lies occur when one or more of <Jp tx or (Jp ty are equal to 
— 1. It appears these anomalies cannot be understood in 
terms of dimensional reduction to d = 1 SPT phases, as 
explained in Sec. V. 

V. DIMENSIONAL REDUCTION VIEWPOINT 

For type 1 and 2 anomalies (see Sec. I), which includes 
the examples described in Sections IV A through IV C, 
the flux-fusion anomaly test can be understood in terms 
of dimensional reduction to d = 1 SPT phases. This 
viewpoint provides a different way of applying the 
anomaly test, which does not depend on some of the 
formalism introduced in Sec. III. In particular, for the 
discussion below, we do not need the description of sym¬ 
metry action on Z„ fluxes f1 presented in Sec. IIIB. Our 
discussion makes significant use of the results of Ref. 25, 
especially for the case of reflection symmetry. 

We imagine putting the un-gauged d = 2 SET phase 
with Z 2 topological order on a cylinder with large but 
finite circumference, in such a way that the symmetry 
is preserved. The longitudinal dimension of the cylin¬ 
der remains infinite. We can approach the limit of an 
infinitely long cylinder either via finite-length cylinders 
with open boundary conditions, or those with their ends 
identified periodically. The minimally entangled states 
(MES) of the SET phase are those with definite anyon 
flux threaded along the cylinder (see Ref. 25 for a more 
complete discussion). The cylinder is a d = 1 system, 
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FIG. 2. Illustration of the flux-fusion anomaly test using di¬ 
mensional reduction to a d = 1 cylinder. The Z 2 C U(l) flux 
17 is threaded twice along the left-hand cylinder, while m is 
threaded along the right-hand cylinder. Due to the fusion rule 
ft 2 = m, these two systems are in the same d = 1 SPT phase. 


gapped and symmetric. Therefore each of the MES is in 
a d = 1 SPT phase. 

Starting from any of the MES, we imagine threading 
an anyon a along the cylinder (see Fig. 2). This can be 
accomplished by a process creating an a-a pair in the 
bulk (a is the anti-particle of a), and then separating a 
and a to infinity. This maps the initial MES to a different 
MES, and so can be thought of as a mapping of d = 1 
SPT phases. In Ref. 25, it was argued that this map¬ 
ping between SPT phases only depends on the nature of 
symmetry action on a in the original d = 2 SET phase. 

Now we consider the case of G = U(l) x Zf symmetry, 
making our usual assumptions that the e-particle (m- 
particle) carries half-odd-integer (integer) U(l) charge. 
Starting with some MES, we thread to along the cylin¬ 
der. We use only the ZJf symmetry to analyze the dimen¬ 
sionally reduced d = 1 SPT phases, so that there are two 
phases distinguished by a Z 2 invariant, that corresponds 
to the presence or absence of Kramers doublet end states 
at the open boundaries [60, 61]. If the to is a Kramers 
singlet [(T m ) 2 = 1], then threading it along the cylinder 
leaves this Z 2 invariant unchanged. On the other hand, 
if to is a Kramers doublet [(T m ) 2 = —1], threading it 
along the cylinder flips the Z 2 invariant. 

To apply the anomaly test via dimensional reduction, 
we next introduce fluxes of the Z 2 C U(l) symmetry. 
Note that we do not fully gauge the Z 2 symmetry; it is 
enough to consider static flux defects of this symmetry, 
without introducing a dynamical gauge field. We can 
thread the flux 17 along the cylinder, which amounts to 
introducing a flux defect 17 near one end, and the cor¬ 
responding anti-defect fl near the other end. Because 
the time reversal symmetry maps 17 1 —> 12 (as shown 
in Appendix B), this can be done while preserving Z^, 
and threading 12 gives another map between d = 1 SPT 
phases. This map either flips the Z 2 invariant or leaves 
it the same. 

Finally, because f l 2 = m , threading f l twice along the 


cylinder is the same as threading to, as illustrated in 
Fig. 2. But threading II twice must leave the Z 2 SPT in¬ 
variant unchanged, which means that threading to must 
also leave this invariant unchanged. This recovers the re¬ 
sult that to being a Kramers doublet is anomalous [62]. 

In this argument, we used the fact that time rever¬ 
sal maps O 1 —>• O. It is important to note that this is 
not a priori obvious simply because 17 is a 7 r flux of the 
U(l) symmetry, because 17 may be attached to an anyon 
under the action of symmetry, as indeed occurs in the 
following example. Instead, we need the considerations 
of Appendix B to conclude 17 1 —> 17 under time reversal. 

For G = U(l) x Zf symmetry, a very similar discus¬ 
sion applies. We choose the Z^ symmetry to exchange 
the two ends of the cylinder; that is, the d = 2 reflec¬ 
tion symmetry becomes d = 1 reflection symmetry upon 
dimensional reduction. SPT phases in d, = 1 protected 
by such Zf symmetry are also characterized by a Z 2 in¬ 
variant [60, 61]. Ref. 25 argued that threading to along 
the cylinder preserves this invariant if (P m ) 2 = l, and 
flips the invariant if (P m ) 2 = — 1 . Again, we consider 
the effect of threading a flux 17, of the Z 2 C U(l) sym¬ 
metry, along the cylinder. Because Zf maps 17 to the 
anti-flux S2 = I7 3 = mI7 (Appendix B), this can be done 
while preserving the Z^ symmetry. Therefore thread¬ 
ing 17 either flips or preserves the Z 2 SPT invariant. At 
this point the discussion proceeds identically to the case 
of time reversal above, and we find that the symmetry 
fractionalization pattern with (P m ) 2 = —1 is anomalous. 

This discussion also applies directly to the case of 
G = U(l) x pm symmetry, because the anomalous sym¬ 
metry fractionalization patterns found in Sec. IV C are 
associated with two different Zff subgroups of pm. One 
of these is generated by P x , and the other is generated 
by T X P X . 

Using the approach of Sec. Ill, we also find anomalies 
associated with the interplay between time reversal and 
other symmetries, that apparently cannot be understood 
from the dimensional reduction point of view. These were 
designated type 3 anomalies in Sec. I, and arise when 
time reversal forms a semi-direct product with U(l) [i.e. 
U(l) xZ[f C G\. They occur in the examples G = (U(l) x 
Z 2 ) x pi (Sec. IVD) and G = (U(l) x Zj) x pm, G = 
(U(l) x Zlf) xjMtoto (Appendix G). For example, in each 
of these cases, T"' (T m ) gives the action of translation in 
the ^-direction (time reversal) on m-particles, and these 
generators obey the relation 

T rn T r x n = a¥ tx T™T m , ( 88 ) 

where we find that <T™ tx = — 1 is anomalous. If we try to 
apply dimensional reduction here, we observe that time 
reversal maps the flux 17 to the anti-flux fl, without ex¬ 
changing the two ends of the cylinder. This means that 
the state obtained upon threading 17 breaks time rever¬ 
sal, and is thus not a d = 1 SPT phase preserving the 
symmetries involved in the anomaly. 
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VI. BOSONIC TOPOLOGICAL CRYSTALLINE 
INSULATORS 

Here, we use the results from the flux-fusion anomaly 
test to identify and distinguish some d = 3 bosonic 
TCIs. We focus on the examples of G = U(l) x Zf 
and G = U(l) xpm symmetry discussed above in Sec. IV; 
the latter example is sufficiently complex to illustrate the 
corresponding general results. The discussion for SPT 
phases with G = U(l) x Zj symmetry entirely parallels 
that given in Sec. VIA, where we thus also comment on 
that case. The focus is on understanding the extent to 
which information obtained from the anomaly test can 
distinguish the TCI phases identified, without using fur¬ 
ther information, while also illustrating what additional 
information can be used to make finer distinctions among 
phases. 

A. G = U(l) x Zf 

In Sec. IV B, we found that the symmetry fractional- 
ization patterns eCmP and eCPmP are anomalous. A 
coupled layer construction (see Ref. 37 and Appendix F) 
shows that these fractionalization patterns - indeed, any 
fractionalization pattern - can be realized as a surface 
of a d = 3 SPT phase, or bosonic TCI. Because each 
fractionalization pattern is anomalous, its corresponding 
SPT phase must be non-trivial. 

Having established that the eCmP and eCPmP SPT 
phases are non-trivial, we would like to understand 
whether these phases are distinct from one another. 
Moreover, as is well-known, SPT phases can be added 
together by combining together two decoupled systems, 
and observing that the combined system thus obtained 
is also an SPT phase. This operation is expected to form 
an Abelian group. We would also like to know how the 
eCmP and eCPmP SPT phases behave under this ad¬ 
dition operation. 

Let us consider adding together two copies of the 
eCmP SPT phase. This results in a surface with two 
decoupled surface SET phases, which we denote by 
eCmP ® eCmP. Denoting the anyons in one SET phase 
by ei,TOi, and in the other by e 2 ,m 2 , we consider the 
result of condensing eie 2 and miTO 2 . This condensation 
destroys the topological order, since all anyons of the 
eCmP © eCmP surface are either condensed, or are con¬ 
fined due to non-trivial mutual statistics with the con¬ 
densate. In addition, both the particles condensed have 
integer charge and P 2 = 1, so that they can be condensed 
without breaking any symmetries. Therefore, we have 
obtained a gapped, symmetry-preserving, trivial surface, 
and the resulting SPT phase is the trivial phase. 

The same conclusion clearly holds for the eCPmP SPT 
phase. Indeed, the conclusion holds for any SPT phase 
with surface Z 2 topological order (at least as long as the 
symmetry does not permute the anyon species). To find 
SPT phases with order higher than two under addition, 


we would need to consider other types of topological or¬ 
der ( e.g. Z n topological order), or potentially those with 
symmetries permuting the anyon species. 

Next, we consider adding together the eCmP and 
eCPmP SPT phases, obtaining a eCmP® eCPmP sur¬ 
face. This surface can be simplified by condensing mim. 2 , 
which again can be done without breaking any symme¬ 
tries. This results in a new surface SET phase with Z 2 
topological order, with anyons e, m, given in terms of the 
original anyons by e = eie 2 and m = mi ~ m 2 . There¬ 
fore, the fractionalization pattern after condensing m\m 2 
is ePmP. 

The flux-fusion anomaly test provides no information 
about ePmP , so without additional information we can¬ 
not draw any further conclusions. It has been shown via 
other methods that ePmP is anomalous [28, 63], so the 
two SPT phases are different. 

The flux-fusion anomaly test on its own allows us to 
distinguish a pair of SPT phases. We can take this pair 
to be either the trivial phase and the eCmP phase, or the 
trivial phase and the eCPmP phase. Both pairs form a 
Z 2 subgroup group under addition of SPT phases; this is 
the result appearing in the fourth column of Table I. 

The same discussion holds for G = U(l) xZj, replacing 
“P” by “T” everywhere, so that we consider the non¬ 
trivial SPT phases with eCmT and eCTmT surfaces. In 
this case, we note that eTmT has also been argued to be 
anomalous [35, 37]. 


B. General results and G = U(l) x pm 

We now consider bosonic TCIs for the general case 
of G = U(l)xG' symmetry, using the example of 
G = U(l) x pm to illustrate the discussion. First, we 
make some fixed but arbitrary choice for the fractional¬ 
ization class of the e particle, [w e ]z 2 £ H 2 (G', Z 2 ), and 
we consider fractionalization patterns of the form 

F l = eC[cj e ]mO[cj % m ], (89) 

where we have introduced an index i to label the distinct 
vison fractionalization classes [w)Jz 2 £ H 2 (G', Z 2 ). Each 
fractionalization pattern F l corresponds to a bosonic TCI 
(d = 3 SPT phase), for which it describes a surface SET 
phase (see Appendix F). In general, not all the F l, s cor¬ 
respond to distinct or non-trivial SPT phases. Formally, 
it will be convenient to refer to a map ip : {F 1 } —> QspTi 
from the set of fractionalization patterns described in 
Eq. (89), to the (Abelian) group of distinct d = 3 SPT 
phases of the given symmetry, which we denote by Qspt- 

Adding together the SPT phases corresponding to Fi 
and Fj gives the surface SET phase Fi ® Fj. Labeling 
the Fi anyons by ei,mi, and the Fj anyons by e 2 ,m 2 , 
this surface theory can be simplified by condensing eie 2 , 
which can be done without breaking symmetry, because 
both e particles have the same fractionalization class. 
The resulting surface theory has Z 2 topological order, 
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with anyons e = e± ~ e 2 and m = mi m 2 , so that we 
have 

Fi © Fj ~ eG[w e ]mO[«]. (90) 

We thus see that the fractionalization patterns add ac¬ 
cording to the same operation governing multiplication 
of [<] Z2 in .ff 2 (G',Z 2 ), and therefore the set {F 1 } can 
be viewed as a group isomorphic to F 2 (G',Z 2 ). In addi¬ 
tion, this shows that the map p can be viewed as a group 
homomorphism tp : H 2 (G', Z 2 ) —► Qspt- 

We know the group H 2 (G', Z 2 ), and the flux-fusion 
anomaly test gives us some knowledge about the map 92 . 
The goal is to use this information to learn as much as 
possible about Qspt- If F l is a fractionalization pattern 
known to be anomalous, then p{F l ) ^ 1; that is, the cor¬ 
responding SPT phase is non-trivial. In the example G = 
U(l) x pm, recall that [w m ]z 2 = (cr% ty , cr™ px , a™ , <r% px ), 
and the corresponding SPT phase is non-trivial whenever 
= - 1 , or a% px = - 1 , or both. 

In addition, we would ideally like to know which frac¬ 
tionalization patterns are non-anomalous, and thus map 
to the trivial SPT phase. Formally, the set of such frac¬ 
tionalization patterns is the kernel of <p, and is a sub¬ 
group of H 2 {G ', Z 2 ) denoted by Ker p. In general, we do 
not know Ker tp. However, we do know which fraction¬ 
alization patterns are anomaly-negative, these also form 
a subgroup denoted TV" C F 2 (G', Z 2 ). In the present 
example, TV" = Z 2 x Z 2 consists of those fractionaliza¬ 
tion classes of the form [w m ]z 2 = (a" xty , cr^ px , 1,1). In 
general, we have Ker p C A/" C F 2 (G',Z 2 ); that is, non- 
anomalous fractionalization patterns are a subgroup of 
anomaly-negative ones. 

Now, we consider the group S = H 2 (G', Z 2 )/7V, ele¬ 
ments of which are cosets of TV". We will see that there 
are at least as many distinct SPT phases as there are el¬ 
ements of S. From the fact Kerp C M C U 2 (G', Z 2 ), it 
follows immediately that distinct elements of S map to 
disjoint sets of SPT phases in Qspt-[ 64] The disjoint sets 
of SPT phases are thus labeled by elements of S; this is 
the group that appears in the fourth column of Table I. 

In the present example, S = Z 2 x Z 2 , and its elements 
Pi,... ,P4 are the four cosets 


Pi = (1,1,1,1) x TV" 

(91) 

P 2 = (1,1, —1,1) x TV" 

(92) 

P3 = (1, 1, 1, —1) x A/" 

(93) 

Pa = (1,1,-1,-1) x A/". 

(94) 


Each of these cosets corresponds to four different surface 
SET phases, depending on the element chosen from TV". 
Surface SET phases belonging to the same coset may or 
may not correspond to the same SPT phase, but surface 
SET phases belonging to different cosets correspond to 
different SPT phases. Therefore, in this example, there 
are at least four bosonic TCIs (one of which is trivial). 

To obtain additional information, we need to deter¬ 
mine Ker p C A f. For example, in some cases it may be 


true that Ker p = A/", if all the anomaly-negative frac¬ 
tionalization patterns are in fact non-anomalous. The 
number of distinct SPT phases obtained from each coset 
of TV" is |TV"|/| Ker 

Throughout this discussion, we have fixed the e par¬ 
ticle fractionalization class, but it is also natural to 
consider fractionalization pattens with different e par¬ 
ticle fractionalization classes, as we did for the case of 
G = U(l) x Zf symmetry in Sec. VIA. Suppose we add 
together F 1 = eC[w)]mO[w)„] and F 2 = eG[w 2 ]? 7 i 0 [wjjJ, 
to obtain a F 1 © F 2 surface, where [wj]z 2 ^ [w 2 ]z 2 . 
If it happens that [ui^] z 2 = [wjlJJz 2 , we can condense 
toito 2 to obtain a surface SET phase with fractionaliza¬ 
tion pattern eOfWgWgJmOfo;^]. Here, none of the anyons 
carry fractional U(l) charge. There is not a general un¬ 
derstanding of which such fractionalization patterns are 
anomalous. However, many such patterns can be explic¬ 
itly constructed strictly in d = 2 , using, for instance, ex¬ 
actly solvable models or parton gauge theory [65]. This 
has been done for square lattice space group symmetry 
using exactly solvable models [65], and could be done for 
other symmetry groups as needed. In addition, in the 
case of reflection symmetry, the ePmP fractionalization 
pattern is anomalous [28, 63]. Results along these lines 
can thus be used to obtain further information on bosonic 
TCI phases, which we leave for future work. 


VII. ANOMALOUS SUPERFLUIDS 

Our results on anomalous symmetry fractionalization 
patterns can also be used to obtain anomalous surface 
superfluid states of d = 3 bosonic TCIs. The anoma¬ 
lous nature of these superfluids arises from the symmetry 
fractionalization of vortices, which transform projectively 
under G' symmetry. It is particularly useful to identify 
such anomalous superfluids, because it is straightforward 
to proceed from their formal description to explicit field 
theories, which can be used to describe not only the sur¬ 
face superfluid phase, but also nearby surface phases and 
phase transitions. 

A related prior work is Ref. 35, which studied d = 3 
bosonic topological insulators (with U(l) and time rever¬ 
sal symmetry) by constructing field theories for anoma¬ 
lous surface superfluids. Some of those superfluids are 
characterized by non-trivial vortex symmetry fractional¬ 
ization, and were argued to be anomalous based on the 
impossibility of realizing a trivial, gapped surface in an 
explicit dual vortex field theory for the surface. Our re¬ 
sults are complementary, allowing one to establish that 
some dual vortex field theories for d = 2 superfluids are 
anomalous without a detailed and potentially subtle anal¬ 
ysis of possible phases. 

As usual, we consider a surface SET phase with Z 2 
topological order and symmetry fractionalization pattern 
eC[u) e ]mO[u) m ], but making the additional assumption 
that [w e ] Z2 = 1 ; that is, the e particle transforms triv¬ 
ially under G'. It is therefore possible to condense an 
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e particle carrying U(l) charge 1/2 and obtain a super- 
fluid, with spontaneously broken U(l) symmetry, where 
G' symmetry is preserved. Under these circumstances, 
the vison of the SET phase becomes the elementary 27T- 
vortex of the superfluid [46], so the vortex thus inherits 
the G' transformation properties of the m particle. If we 
start with an anomalous surface SET phase, the resulting 
surface superfluid must also be anomalous, because both 
are surfaces of the same non-trivial SPT phase. 

It is well-known that vortices can transform projec- 
tively under symmetry [53-55]. This can be seen con¬ 
veniently in the dual description of a superfluid, where 
vortices carry the charge of a non-compact U(l) gauge 
field, for which the photon is nothing but the superfluid 
sound mode. Symmetry operations acting on vortices can 
thus be augmented by U(l) gauge transformations, and 
the symmetry acts projectively. In fact, in more detail, 
vortices transform as a f-twisted U(l) projective repre¬ 
sentation of G' , and vortex fractionalization classes are 
thus elements of H 2 (G\ U(l)). This can be seen by in¬ 
troducing field operators for the vortices, as was done in 
Sec. IIIB for the Z n flux f2. Here, the field operators 
are labeled by an integer, which is simply the vortic- 
ity, and we have U(l) gauge transformations rather than 
Z 2 n gauge transformations. Otherwise, the discussion 
entirely parallels that given in Sec. IIIB. 

Because a m particle becomes a vortex upon entering 
the superfluid phase, the vortex fractionalization class 
[w«] G Hf(G', U(l)) is given by 

[w„] = [w m ] U( i) = P2([u>m] Z 2 ). (95) 

Remarkably, [w m ]u(i)> which provides a simple mathe¬ 
matical characterization of which fractionalization pat¬ 
terns are anomaly-negative, also has direct physical 
meaning as the fractionalization class of vortices in the 
superfluid phase. This allows us to establish that su¬ 
perfluids with certain vortex fractionalization classes are 
anomalous. 

This conclusion is bolstered by proceeding in the re¬ 
verse direction; that is, we can start with a superfluid, 
and condense pairs of vortices to obtain a SET phase with 
7*2 topological order. This can be done without breaking 
symmetry as long as the vortex fractionalization class 
satisfies [w „] 2 = 1 , so that vortex pairs transform triv¬ 
ially. The fractionalization class of the m particle must 
satisfy Eq. (95), but we note this does not completely 
determine [w m ]z 2 given [w„]. We expect that, given [w„], 
the different possible choices of [w m ]z 2 correspond to in¬ 
equivalent condensates of paired vortices; detailed study 
of this point is left for future work. 

Which vortex fractionalization classes are anomalous? 
We answer this question for the example of G = U(l) x 
pm symmetry, and then make some comments on the 
answer more generally. First, Eq. (95) implies that [uj v ] = 
{atxty = ±1 ,a px ,atxpx) is anomalous whenever a px = 
— 1 , ottxpx = — 1 , or both, because these [u> v \ are obtained 
from anomalous [w m ]z 2 - 


We can also find more anomalous vortex fractionaliza¬ 
tion classes, by starting with an anomalous superfluid, 
and adding a layer of d = 2 superfluid. This can be done 
by first assuming that each layer has an independent U(l) 
symmetry, and then breaking the resulting U(l) x U(l) 
down to U(l); that is, we allow unit charge excitations 
to tunnel between the two layers. Before breaking the 
U(l) x U(l) symmetry, each layer has independent vor¬ 
tices, schematically labeled by iq and u 2 . After breaking 
the symmetry, vi and V 2 vortices are confined together, 
so that the new superfluid state has vortices v = v±V 2 - 
This results is a modified vortex fractionalization class 
[<jj v ] = [x’ 7 ,| ] [uj V2 ]. 

In the present example, vortex fractionalization classes 
[w„] = (a txty , 1,1), with a txty an arbitrary U(l) phase, 
can occur in d = 2. Writing a tx t y = e 2nm , such su¬ 
perfluids occur for bosons on the square lattice at filling 
n. This is easily seen via straightforward application of 
boson-vortex duality to such a model; briefly, the vor¬ 
tices feel the background boson density as a magnetic 
flux of 2irn per plaquette, and thus transform projec¬ 
tively under translation symmetry. By adding layers 
of such non-anomalous superfluids, we can see that the 
only non-anomalous vortex fractionalization classes are 
[w„] = ( atxty , 1 , 1 ), and all others are anomalous. 

This result can be stated in a more general fash¬ 
ion, namely, [oj t ,] is non-anomalous if and only if it 
can be continuously deformed to the identity element of 
H 2 (G', U(l)). We conjecture that this result holds for all 
symmetries G = U(l) x G', but we do not have a general 
argument, for two reasons. First, we note that the discus¬ 
sion above relied on being able to find all non-anomalous 
vortex fractionalization classes for G = U(l) x pm. sym¬ 
metry. Second, in this case, H 2 (G', U(l)) was a prod¬ 
uct of U(l) and Z 2 factors, containing no Z n factors for 
n > 2. If there were such Z„ factors, some vortex frac¬ 
tionalization classes could not be obtained by condensing 
the e particle of a SET phase with Z 2 topological order. 
It is likely that this could be handled by generalizing the 
flux-fusion approach to SET phases with Z„ topological 
order, a problem that is left for future work. 

We conclude this section, and illustrate the utility of 
the present results, by describing the construction of a 
field theory for the surface of a bosonic TCI with sym¬ 
metry G = U(l) x pm. We work in a dual description 
of the surface superfluid, introducing a two-component 
complex field <f>.„ for the superfluid vortices. carries 
unit charge of the U(l) gauge field a M , in terms of which 
the global U(l) current is = e pv \d v a\/2 tt. We work 
in Euclidean space time. We choose the pm symmetry 
generators to act on the vortices as follows: 

T x ,T y : $ v (x,y,r) $ v (x,y,T) (96) 

P x : $ v (x,y,r) -> (ia y )^ v (-x,y,r). (97) 

As usual, we neglect the action of lattice translations T x 
and T y on the spatial position of the continuum field , 
as this only leads to subleading gradient terms. The pres¬ 
ence of the Pauli matrix ia y in the action of P x implies 
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that = — 1 acting on $ t ,, and we have the vortex 
symmetry fractionalization [w„] = (1, —1, —1). This frac- 
tionalization class is anomalous, so we are not describing 
a d = 2 superfluid, but rather the surface of a bosonic 
TCI. We note that we chose as a two-component field 
in order to realize this non-trivial fractionalization class. 

The Lagrangian is 

C=\{dn+ ia^y.b v \ 2 + + U(4>„) H-. (98) 

A* 

Here, the first term is the vortex kinetic energy, the sec¬ 
ond term controls the superfluid stiffness, and V($„) is 
a potential for the vortex field, whose form is dictated by 
gauge invariance and the action of the microscopic sym¬ 
metries. The ellipsis includes various other perturbations 
allowed by symmetry. This field theory can be used to 
study the superfluid phase itself (where d>.„ is massive), 
neighboring phases described as condensates of 4> t , (which 
break lattice symmetries), surface SET phases where 
is condensed, and transitions among these phases. 

VIII. DISCUSSION AND OUTLOOK 

We introduced the flux-fusion anomaly test, a method 
to detect some anomalous symmetry fractionalization 
patterns in d = 2 SET phases with Z 2 topological or¬ 
der. This constrains the possible physical properties of 
strictly d = 2 SET phases, and is a step toward the 
full classification of such phases in the presence of crys¬ 
talline symmetry. In addition, the same results allow us 
to identify and distinguish some d = 3 SPT phases via 
their surface theories, including some bosonic TCIs that 
have not previously been identified to our knowledge. For 
some of the bosonic TCIs, we identified not only surface 
SET phases with anomalous symmetry fractionalization, 
but also anomalous surface superfluids distinguished by 
the projective symmetry transformations of vortex exci¬ 
tations. 

We note that the flux-fusion anomaly test is closely re¬ 
lated but not equivalent to a “monopole tunneling” ap¬ 
proach developed in [36] and used in [37]. In this ap¬ 
proach, one considers a d = 3 SPT phase with U(l) sym¬ 
metry, gauges the U(l) symmetry, and studies magnetic 
monopole excitations in the bulk. If one has U(l) x Z^ 
symmetry, the monopole can be a Kramers doublet, in¬ 
dicating the bulk SPT phase is non-trivial. Tunneling 
a monopole into the bulk through a superfluid surface 
leaves a vortex behind on the surface, which must also 
be a Kramers doublet. Condensing double vortices on the 
surface leads to the eCmT state, which is thus a surface 
SET of a non-trivial SPT phase. 

Using the description of symmetry action on vortices 
given in Sec. VII, very similar reasoning can be applied 
to the symmetries considered in this paper, and the same 
anomalies we find can presumably be diagnosed. How¬ 
ever, the two approaches are not equivalent. In particu¬ 
lar, the flux-fusion method is more general as it can be 


applied for discrete symmetries, e.g. G = Z n x G 7 , as 
mentioned below. 

In this paper, we focused primarily on symmetries of 
the form G = U(l) x G space and G = [U(l) x T\] x G S p ace , 
where G space is a d = 2 space group. The latter symme¬ 
try is particularly relevant for systems of bosons. We 
did not consider the very important class of symmetries 
G = U(l) x Zjf x G space , which are relevant for spin sys¬ 
tems with continuous spin rotation symmetry. For exam¬ 
ple, when U(l) C SO(3), these are the symmetry groups 
of Heisenberg spin models. The reason for this omis¬ 
sion is a surprising finding that complicates application 
of our anomaly test: for these symmetries, it is some¬ 
times impossible to gauge Z n C U(l) without breaking 
Z^ x G spaC e symmetry, even in strictly two dimensions 
[40]. This can occur when some lattice sites transform in 
a projective representation of the on-site U(l) x Zif C G 
symmetry; for example, when U(l) C SO(3), this means 
that there are S = 1/2 or other half-odd-integer spins in 
the system. Naively, it would appear the anomaly test 
is less useful for these symmetries, but, remarkably, it 
turns out that this obstruction to gauging the symme¬ 
try makes the anomaly test significantly more powerful. 
These results will be presented in a forthcoming work 

[40]. 

More generally, to which symmetry groups does the 
flux-fusion anomaly test apply? One point is that U(l) 
symmetry is not required, and it is simple to generalize 
the results of this paper to symmetries G = Z n x G' in a 
straightforward manner. This works as long as n is even, 
so that we can sensibly choose e to carry half-charge of 
Z n , and as long as the G 7 symmetry constrains the Z„ 
flux to be a boson. We note that, if G = Z 2 x G 7 , the Z 2 
flux is always a boson, independent of G 7 (see Sec. Ill A). 
In principle, we can also consider G = G 0 x G 7 , where 
G 0 is some finite, on-site, unitary symmetry. In practice, 
in the latter case, one generally obtains a non-Abelian 
gauged SET phase, which can be expected to increase 
the complexity of analysis required. 

A related point is that the requirement that symmetry 
fluxes are bosons is not fundamental, but is rather im¬ 
posed because it simplifies the analysis. For on-site, uni¬ 
tary symmetries that do not permute the anyons of the 
gauged SET phase, we believe it likely this requirement 
plays no role and can simply be ignored. More generally, 
one needs a description of the action of symmetry on non- 
bosonic anyons, which is subtle and not yet fully under¬ 
stood for crystalline symmetry [19, 21, 25, 28]. However, 
we expect that the necessary theory will become available 
with further progress, in which case it can be applied to 
broaden the applicability of the flux-fusion anomaly test. 

It is also interesting to consider generalizing the flux- 
fusion anomaly test to other topological orders. The ba¬ 
sic idea of the anomaly test is, given an action of sym¬ 
metry on the anyons of an un-gauged SET phase, to 
determine whether this action can be extended consis¬ 
tently to symmetry fluxes. This idea applies more gen¬ 
erally to SET phases with topological orders and sym- 


18 


metries beyond those considered here, although we do 
not expect our detailed analysis to apply in general. For 
on-site, discrete, unitary symmetries, the framework of 
G-crossed tensor category provides a comprehensive de¬ 
scription of SET phases [32] and a systematic means of 
detecting anomalies [32, 38]. For symmetries where both 
approaches apply, the flux-fusion anomaly test as devel¬ 
oped in this paper is certainly less general than the G- 
crossed tensor category approach, but it has the advan¬ 
tage of identifying some anomalies in a physically intu¬ 
itive way. Moreover, without the need to introduce fluxes 
for all symmetries, the flux-fusion approach can be easily 
applied to continuous, anti-unitary and spatial symme¬ 
try, as illustrated by the examples discussed in this paper. 

The examples studied in this paper can be analyzed 
without resort to G-crossed tensor category theory due 
to the simplicity of the topological orders involved. The 
magnetic sectors of the Z 2 topological order and the 
gauged Z 2n theory have trivial F and R matrices. There¬ 
fore, when analyzing their transformation under symme¬ 
try, we do not need to worry about the “gauge trans¬ 
formation” on fusion spaces, as discussed in Eq. 58 of 
Ref. 32. This greatly simplifies the mathematical struc¬ 
ture involved, and the flux fusion procedure as discussed 
in this paper can be implemented. 

General SET phases can have nontrivial F and R ma¬ 
trices, and it is important to take “gauge transforma¬ 
tions” into account when analyzing symmetry action. To 
avoid this complexity, we can restrict to the case where 
the symmetry flux f lh (for h £ G) remains invariant un¬ 
der the symmetry action of g £ G. That is, we require (1) 
h commutes with g , so that remains the same symme¬ 
try flux, and also (2) Qh remains in the same topological 
sector and is not attached to an anyon under the ac¬ 
tion of g. The second condition can be violated when 
g and h act non-commutatively on the anyons. For ex¬ 
ample, in the projective semion example of Ref. 38, with 
G = Z 2 x Z 2 symmetry, the two Z 2 symmetries anti¬ 
commute with each other on the semion, and the flux 
of one Z 2 is glued to a semion under the action of the 
other Z 2 . The F and R matrices are non-trivial in this 
example, so we do not expect a straightforward general¬ 
ization of the flux fusion method described in this paper 
to apply. 

When the above two conditions are satisfied, symmetry 
g has a local action on and we can talk about the 
symmetry fractionalization of g on f R without worrying 
about “gauge transformations.” Here we remark that 
g and h can be the same type of symmetry operation. 
More precisely, in the main text we only discussed cases 
where g and h lie in two different factors of a semidirect 
product. However, this is not necessary, and the flux 
fusion idea can apply even when g = h. We discuss such 
an example in a study of d = 3 SET phases [66] . 

Beyond the anomaly test itself, one natural direction 
for further studies is to develop an understanding of the 
physical properties of the bosonic TCIs we have identi¬ 
fied. In light of prior work on bosonic topological insu¬ 


lators with U(l) and time reversal symmetry [35], we 
expect that the surface dual vortex field theories dis¬ 
cussed in Sec. VII will be particularly useful in this re¬ 
gard. Along the same lines, it will be interesting to look 
for simple, physically reasonable models realizing bosonic 
TCIs. 

We also hope that our results on bosonic TCIs will be 
useful as a stepping stone to identify and perhaps classify 
TCIs of interacting electrons. As has been established 
for electronic topological insulators [with U(l) and time 
reversal symmetry], there are non-trivial electronic topo¬ 
logical phases that can be understood by forming com¬ 
posite bosonic particles out of electrons (Cooper pairs, 
or spins), and putting these objects into a bosonic SPT 
phase [67]. This is an important part of the classifica¬ 
tion of interacting electronic topological insulators given 
in Ref. 67. 

Note added. During the review process of this pa¬ 
per, some closely related work has appeared. In par¬ 
ticular, Ref. 68 extended the flux-fusion anomaly test 
to Z 2 spin liquids with SO(3) spin rotation symmetry, 
and showed that the vison symmetry fractionalization in 
5=1/2 Heisenberg models on square and kagome lat¬ 
tices is completely fixed. References 69 and 70 adapted 
and used flux-fusion to constrain the symmetry fraction¬ 
alization of the chiral spin liquid phase of the kagome 
Heisenberg model. Another related development is the 
work of Ref. 63, which presented an approach to classify 
SPT phases protected by point group symmetry based on 
a kind of dimensional reduction; we anticipate this ap¬ 
proach can be generalized to provide an alternate charac¬ 
terization and more complete classification of the bosonic 
TCIs identified here. 
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Appendix A: Procedure to gauge 7L n c u(i) 
symmetry 

Here, we consider symmetry groups G = (U(l) x Z^) x 
G s , where G s is a d = 2 space group, and describe an 
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explicit procedure to gauge the Z n C U(l) symmetry. In 
particular, we verify that this can be done while preserv¬ 
ing G' = Z^ x G s symmetry. We also discuss the case of 
G = U(l) x Zj, giving a procedure to gauge Z n C U(l) 
while preserving Z|\ While these conclusions may appear 
obvious, they do not hold in general for other symmetry 
groups (in particular, for G = U(l) x Z% xG s ). This has 
interesting consequences that will be explored in a future 
publication [40]. 

First we discuss the case G = (U(l) x Z^) x G s . We 
consider a bosonic model defined on a lattice with sites 
r, which is invariant under the space group symmetry 
G s . Each g £ G s acts on lattice sites, which we write 
formally as n-> gr. There is a Hilbert space associ¬ 
ated with each lattice site, and the full Hilbert space is 
the tensor product T~L = (g ) r T~L r ■ Because the U(l) sym¬ 
metry is on-site, for each lattice site r there is a charge 
density operator N r with integer eigenvalues. Because 
time reversal forms a semidirect product with the U(l), 
we have 


TNrT- 1 = N r . (Al) 

In general, we might wish to allow for a shift N r —> 
N r + SN r under time reversal. But, since we assume 
the ground state is invariant under T, we must have 
(N r ) = ( N r ) + 5N r , and 8N r = 0. Moreover, the space 
group operation g £ G s is represented by U g , and acts 
on the charge density by 

UgNrU- 1 = N gr . (A2) 

To gauge the Z n C U(l) symmetry, we introduce Z n 
electric field and vector potential operators, that reside 
on oriented links £ = (r, r'), where each link joins a pair 
of lattice sites r and r'. The set of links is chosen to 
make the lattice into a connected graph that is invariant 
under space group symmetry; for example, choosing links 
to join nearest-neighbor sites is sufficient in many cases. 

The electric field et and vector potential an act on the 
Hilbert space of link £ = (r,r'), which we choose to be 
n-dimensional with basis {10), |1),..., \n — 1)}. The link 
operators are defined by 


a e \k) - exp ( ^ \k) 

(A3) 

ei\k) = (fc + 1) mod n). 

(A4) 

These lattice vector fields are oriented, so that if £ = 
(r',r) is £ with reversed orientation, then ei = ej and 

ag = a\. 

We impose the Gauss’ law constraint 


-j—r [2ni 

\\ee = exp - N r , 

ir^r 

(AS) 

where the product is over those links £ that join 
sites, with orientation pointing away from r. 

r to other 
Choosing 


time reversal and space group operations g £ G s to act 
on e rr / by 

Terror- 1 = e\ r , (A6) 

UgC rr 'Ug = Cgr,gr' •> (AT) 

we see that the Gauss’ law constraint respects the G' 
symmetry. In addition, the Hamiltonian has to be made 
gauge-invariant via the minimal coupling prescription, 
which can be done while respecting G'. 

Now we consider the case G = U(l) x Z|\ Again we 
have a lattice with sites r; because there is no space group 
symmetry, the lattice does not have to satisfy any sym¬ 
metry conditions. Each site again has a charge density 
operator N r with integer eigenvalues. Time reversal now 
acts by 

TNrT- 1 = 8N r - Nr, (A8) 

where 5N r must be an integer. By making constant in¬ 
teger shifts of N r , we may choose 5N r = 0,1. Next, by 
combining pairs of sites together as needed, and making 
further integer shifts of N r , we can set SN r = 0. 

Introducing Z n electric fields and vector potentials as 
above, and imposing the Gauss’ law Eq. (A5), we choose 
time reversal to act on the electric field by 

Te rr 'T _1 = e rr '- (A9) 

With this choice, the Gauss’ law constraint respects the 
Z^ symmetry, as desired. 


Appendix B: Conditions under which 9 is a boson, 
and permutation of anyons in the gauged SET phase 


Here, we show that the Z n symmetry flux H is a bo¬ 
son in the gauged SET phase whenever time reversal or 
reflection symmetry is present. We also show that these 
operations either map H H>■ H, or H i —> H 2n_1 , depending 
on whether they commute with the U(l) symmetry. 

The starting point for the analyses below are the fusion 
rules and statistics of the gauged SET phase. According 
to the discussion of Sec. Ill A, the fusion rules are 


Q n = 1 
e 2 = Q 
m 2 = 1 

n n = mQ k , 
and the statistics are specified by 
0 e = 0 m = 0 
= 7T 

9 q = ©e,Q = ©m,Q = 0 

2-7T 


©Q ,0 = 


n 
7 r 


@e,n = 

n 

0m, n = 0 

o« = T 

n z 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 

(B7) 

(B8) 

(B9) 

(B10) 

(BH) 
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Here, 0 < k < n is an even integer. 

The statistics of the gauged SET phase must obey cer¬ 
tain conditions in the presence of time reversal or reflec¬ 
tion symmetry. For Zj time reversal symmetry gener¬ 
ated by T, we write the action of T on some anyon a 
in the gauged SET phase as T * a. The statistics must 
satisfy 

e r *a = -e a (B12) 

0T*a,T*b = — 0a,b- (B13) 

These relations hold because the time reversed (clock¬ 
wise) exchange process with time reversed anyons must 
give the same result as the ordinary exchange process 
before time reversal. 

Next, under Zf reflection symmetry generated by P, 
we denote the action of P on a by P*a. Equations (B12) 
and (B13) again hold, simply replacing T by P. This is 
the case because a counterclockwise exchange process is 
mapped to a clockwise one under P. 

We will use these relations to show that H is a boson 
whenever Z^ or Z^ symmetry is present. There are four 
cases to consider, where G contains a subgroup U(l) x 
Z%, U(l) x Zf\ U(l) x Tl or U(l) x Zf. We handle these 
cases one by one: 

Case 1 : G contains a subgroup U(l) x Z|\ 

Because T reverses the sign of U(l) charge, in the 
gauged SET phase we have 

T*Q = Q = Q n ~ l . (B14) 

To determine the action of T on e and to, note that 
T leaves these anyons invariant in the un-gauged SET 
phase, but it reverses their U(l) symmetry charges. The 
e sector of the gauged SET phase consists of those e par¬ 
ticles of the un-gauged SET phase whose U(l) charge 
modulo n is 1/2. Similarly, the to sector in the gauged 
SET phase consists of those to particles in the un-gauged 


This implies q = 1, and so far we have shown Cl' = QPQ. 

Finally, we consider the self-statistics of Q, and apply 
Eq. (B12), finding 


^=0 n = -0 a . = - 0 , 

= -p®Q,n 
irk 2irp 
n 2 n 


Q?n 


This implies 


27 rk 
n 2 


2irp 

1 

n 


which has the unique solution k = p = 0. 
Therefore we have shown that 

r*ci = ci. 


(B22) 

(B23) 

(B24) 

(B25) 

(B26) 


We also showed that k = 0, so that His a boson (6q = 0). 
Case 2: G contains a subgroup U(l) x Zf. 

In this case, reflection does not act on U(l) charge, 
and we have 


P*Q = Q 

(B27) 

P * e = e 

(B28) 

P * TO = TO. 

(B29) 

T * H and write 

H' = eW, 

(B30) 


for integers 0 < p 1 q < 2n — 1 to be determined. 

We follow the same strategy as in Case 1, repeatedly 
applying Eqs. (B13) and (B12) to determine Cl'. First we 
have 


SET phase with the U(l) charge Omodn 
have 

. Therefore, we 

0 — — ®m.Q' — i 

(B31) 

T* e = Qe 

(B15) 

which implies p is even. We define p = p/ 2, 

so that 

T * TO = TO. 

(B16) 

W = QPW. Then 


Now, let Cl! = T * H. In general, we can write 

7 r 

— = ©e,Q = —O e ,Q' 

(B32) 

H' = eW, 

(B17) 

n 

7 r 


for integers 0 < p, q < 2n — 1 that we 

will determine. 

= ~q®e,n = ~q~- 

n 

(B33) 


This is a unique parametrization of all (2?i) 2 anyons in 
the gauged SET phase. 

Using Eq. (B13), 

0 = 0m,n = 0771,0' = P®m,e = p7T. (B18) 

This implies p is even, and letting p = p/2, we have 
il' = QPtt q . Next, we apply Eq. (B13) again, this time 
to the mutual statistics of e and H, to obtain 


This implies that q = 2n — 1, and so far we have shown 
H' = QPQ 211 - 1 . Finally, 


^ = o n = -e n > 


= ~ [^n 2 ™- 1 + ©qp.o 2 "- 1 ] 


(B34) 

(B35) 


- = 0 e ,n = -0 
n 


eQ, O' — ~®eQ n ~ 1 ,Qp O^ 


TT „ 27T 

q- + {n- 1 )q — 

. n n J 


7T 

= 5 —• 
n 


(B19) 

(B20) 

(B21) 


= -[(2n - 1) 2 ^| + (2n- l)p —]. (B36) 

rr n 

Rearranging terms, and dropping those that vanish mod¬ 
ulo 27 t, this is equivalent to 


2nk 2 ,tt 
— = — (p + 2k), 


(B37) 
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which has the unique solution k = p = 0 . 

Therefore we have shown that 

P*n = n 2n ~ 1 . (B38) 

We also showed that k = 0, so that 12 is a boson (On = 0). 
Case 3: G contains a subgroup U(l) x Z|\ 

Here, time reversal does not change the U(l) charge, 
so we have 

T*Q = Q (B39) 

T * e = e (B40) 

T * m = m. (B41) 

These equations are identical to those for P in case 2. Be¬ 
cause the symmetry conditions on statistics are the same 
for T and P symmetry, the analysis proceeds exactly as 
in case 2, and we have T* 12 = 12 2ra_1 and On = 0 . 

Case 4 : G contains a subgroup U(l) x Zf. 

In this case, P reverses U(l) charge, so as in case 1 we 
have 


V*Q = Q 

(B42) 

V -ke = Qe 

(B43) 

V -km = m. 

(B44) 


Because these equations are identical to those in case 1, 
the analysis proceeds identically, so that P * 12 = 12 and 
On = 0. 


Appendix C: Specifying fractionalization classes in 
terms of U(l) and G' fractionalization classes 

By definition, the fractionalization class of e or m is 
an element [w] £ H 2 (G, Z 2 ). In this paper, we consider 
G = U(l)xG', and we specify the fractionalization class 
by two pieces of information: 1 ) whether the particle 
carries integer or half-odd integer U(l) charge, and 2) 
an element [<*/] £ H 2 (G\ Z 2 ). Here, we show that all 
fractionalization classes can be uniquely specified in this 
manner. 

We observe that [w] £ H 2 {G, Z 2 ) uniquely deter¬ 
mines elements of P 2 (U(1),Z 2 ) [corresponding to the 
U(l) charge modulo 1] and H 2 (G', Z 2 ). These elements 
are obtained by restricting the arguments of the factor 
set oo(gi,g 2 ) to the U(l) and G' subgroups, respectively. 
Therefore, we need only show that no additional infor¬ 
mation is needed to uniquely specify [w]. 

We consider a projective representation of G , where 
4> £ U(l) is represented by e 1 ^ , and g £ G' is repre¬ 
sented by T(g). Letting a q = 1 (a q = —1) correspond to 
integer (half-odd-integer) U(l) charge, we have 

e 2nlQ = cr q (Cl) 

r( 5 i)T(g 2 ) = u)\g 1 ,g 2 )r(g 1 g 2 ). (C2) 

So far, we have only specified a q and [a/] £ H 2 (G', Z 2 ). 

To complete the description of the projective repre¬ 
sentation, we need to describe the multiplication of an 


element of U(l) with an element of G'. First, fix g £ G', 
and suppose that <pg = gcj) for all <f> £ U(l). Then, in the 
projective representation 

e i</>Q r{g)e -i*Q = f g ^r(g), (C3) 

where / g (^>) £ {±1}- Setting 0 = 0, clearly f g ( 0) = 1. 
Moreover, the left-hand side is a continuous function of 
</>, so f g (4>) must also be continuous, and f g (4>) = 1 for 
all (j). 

The other possibility we need to consider is a fixed 
g £ G' where (j>g = g (—</>) for all $ £ U(l). In the 
projective representation, 

e foQ r(ff)e foQ = f g ((j))T{g). (C4) 

Here, the same arguments show that f g ((f>) = 1. 

We have thus shown that the fractionalization class 
[w] £ H 2 {G, Z 2 ) is completely specified by a q and [u/] £ 
H 2 (G', Z 2 ). 


Appendix D: Characterization of anomaly-negative 
fractionalization patterns 

We recall that, by definition, the symmetry fractional¬ 
ization pattern eC[w e ]»n 0 [u m ] is anomaly-negative if and 
only if, for each even n > 2, there exists a 2-twisted Z„ 
factor set 4 l n(gi,g2) solving the equation 

u m {gi, 9i) = [(i>n{gi, 92)] n , (Dl) 

where gi,g 2 £ G'. In this section, we prove a simple 
characterization, stated as Theorem 1 in Sec. IIIB, of 
which m particle fractionalization classes [w m ]z 2 give rise 
to anomaly-negative fractionalization patterns. 

It will be convenient to relate the Z 2 and Z 2n fac¬ 
tor sets Um and </>„ to 2 -twisted U(l) factor sets. To 
proceed, let k be a positive integer. Z 2 (G',Z 2 k) is the 
Abelian group of f-twisted Z 2fc 2-cocycles (factor sets) 
for the group G 1 . B 2 {G' ,Z 2 ^) is the corresponding 

Abelian group of 2-coboundaries, which are factor sets 
of the form uj(gi, g 2 ) = A(gi)[A(g 2 )] t(si) [A( 5 ig 2 )] _1 , for 
A (g) £ Z 2 k- The second cohomology group is defined by 
H 2 (G\ Z 2k ) = Z 2 (G',Z 2 k )/B 2 (G',Z 2k ). There is a pro¬ 
jection homomorphism ir 2 k : Z 2 (G\ Z 2k ) —> H 2 (G', Z 2k ). 
Note that if k = 1, we can drop the t subscripts ev¬ 
erywhere, since in that case the twisting by t(g) is triv¬ 
ial. The same definitions hold for U(l) coefficients, in 
which case we call the projection homomorphism 7 Tu(i) : 
Z 2 (G\ U(l)) -A H 2 (G\ U(l)). There is an obvious in¬ 
clusion map i 2k : Z 2 (G', Z 2k ) -A Zf(G', U(l)), which 
just expresses the fact that oj £ Z 2 (G\ Z 2k ) can also 
be viewed as a t -twisted U(l) factor set. 

For each fc, we would like to define a homomorphism 
p 2k : H 2 (G' ,Z 2k ) -A H 2 (G', U(l)), so that the following 
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diagram is commutative: 

Z 2 t (G',Z 2k ) Z 2 (G', U(l)) 

| 7T 2 fc I ' (D2) 

H 2 (G',Z 2k ) H 2 (G', U(l)) 

In fact, we will see that p 2k is the unique homomorphism 
making this diagram commutative. 

Why do we want to define p 2k l Given u> £ Z 2 (G' ,Z 2k ), 
we can define a U(l) fractionalization class by [oj]xj(i) = 
7 Tu(i)(* 2 fc(w)) £ Hf(G', U(l)). If we can find a unique 
p 2k , commutativity of the diagram tells us that [w]u(i) 
depends, in a unique way, only on the Z 2k fractionaliza¬ 
tion class [uj]z 2k = 7 r 2 k (co) £ Hf(G', Z 2k ), by [w] u(1) = 
P2k{[u]z, 2k ) ■ Therefore, it is meaningful to talk about 
Mu(i) as a function of [ uj]j, 2k ■ 

We define p 2k as follows. Pick some element c £ 
H 2 (G' ,Z 2k ). Choose a representative uj £ Zf(G',Z 2k ) 
so that c = ir 2k (ui). Then define p 2 k {c) = 7 Tu(i)(* 2 fc(w)). 

First, we have to check p 2k is well-defined, which means 
it must be independent of the particular representative ui. 
It is easy to see that w, oJ £ Z 2 (G', Z 2k ) belonging to the 
same cohomology class, also belong to the same cohomol¬ 
ogy class after mapping under i 2k to Z 2 (C, U(l)), so p 2k 
is well-defined. Next, we check p 2k is a homomorphism. 
Let c, d £ H?(G\ Z 2k ), with corresponding representa¬ 
tives We have cd = TT 2k (co)n 2k (ijj') = Tr 2k (uju>'), so 

uiui' is a representative for cd. Then 

P2k(cc) = 7 r u( 1 )('i 2 fc(wu/)) = p 2k (c)p 2k (d). (D3) 

Finally, we check p 2k is the unique homomorphism 
making the diagram commutative. Suppose p 2k also 
makes the diagram commutative, but for some c £ 
Hf(G', Z 2fc ), we have p 2 k (c) ± p 2 k(d). Let w £ 
Z 2 (G' ,Z 2k ) be a representative for c, then it follows that 
7 Tu(i)(* 2 feM) = P 2 fe( 7 T 2 fe(w)) = p 2k (^2 k M), which im¬ 
plies p 2k (c) = p2k(c), a contradiction. 

Now we return to Eq. (Dl). Viewing uj m and (f> n as 
U(l) factor sets, and applying 7Tu(i) to both sides of the 
equation, we have 

l w m]u(l) = ([0n]u(l) ) n ■ (D4) 

Therefore, another way of putting the anomaly test is 
that, in order for [w m ]z 2 to be anomaly-negative, [w m ]u(i) 
must have a nth root in ff 2 (G',U(l)) for all even n > 2. 

Now we can prove the desired characterization 
of anomaly-negative [u) m \ z 2 - We assume that 
Lf 2 (G',U(l)) = U(l) fc x A, where A is a finite product 
of finite cyclic factors. This assumption is true in all the 
examples we studied, and we believe it is likely to hold 
in general. 

Proposition 1 . Suppose that [w m ]u(i) lies in the con¬ 
nected component of Hf(G ', U(l)) that contains the iden¬ 
tity element. Suppose also that H 2 (G', U(l)) = U(l) fc x 
A, where A is a finite product of finite cyclic factors. 
Then \ui m ]z 2 is anomaly-negative. 


Proof. It follows from the assumptions that [w m ]u(i) has 
an nth root in 7L 2 (G', U(l)) for any n > 0. This holds by 
the assumption that H^(G', U(l)) is a product of U(l) 
and finite cyclic factors, so that the connected component 
containing 1 is just a product of U(l)’s. Then we have 

Wm( 5 i, 3 2 ) = A ~ 1 (gi)[X(g2)}~ t{ 9 l) X{g 1 g 2 ){n n (gi 1 g2)) n , 

(D5) 

where £ Z 2 (G',U(1)), and A (g) £ U(l). 

We choose 0 < 6 (g) < 2 ir so that A (g) = e ie ^ 9 \ Then 
we define a(g ) = and we choose 

<t>n(9\,92) = 0!~ 1 (gi)[a(g 2 )}~ tl ' 9 l) a(g 1 g 2 )n n (g 1 ,g 2 ). 

(D 6 ) 

This is by construction a nth root of u> m (gi, ff 2 ), so 
*pn(gi192) £ Zj 2n . To show f> n £ Zf(G\ Z 2n ), we note 
that <f> n (gi,92) clearly satisfies the relevant associativity 
condition. For the given u> m , we have thus constructed a 
solution to Eq. (56) for each even n > 0. □ 

The converse of Proposition 1 is also true: 

Proposition 2 . If [w m ]z 2 is anomaly-negative, and if 
Hf(G’, U(l)) = U(l) fc x A, where A is a finite product 
of finite cyclic factors, then [w m ]u(i) lies in the same 
connected component of Hf(G', U(l)) that contains the 
identity element. 

Proof. Under the assumption, Eq. (D4) holds; that is, 
[w m ]u(i) has an nth root in iL 2 (G',U(l)) for all even 
n > 0 . 

We write [w m ]u(i) = (a, (3), where a £ U(l) fe and f3 £ 
A. We will show that j3 = 1, which implies [w m ]u(i) lies 
in the connected component containing the identity. 

Write A = Z P1 x ••• x Z PN , and /3 = (61,... , bjv). 
Observe that [u m ]L = 1, which implies /3 2 = 1. If pi is 
odd, j3 2 = 1 implies 6 ,; = 1. Now consider pi even. Then, 
by assumption, there exists a pA h root of j3, /3 = 7 Pi for 
7 £ A. This implies bi = c Pi for some c £ Z Pi , but for 
any c £ Z Pi , bi = c Pi = 1. Therefore, /3 = 1. □ 

Taking these two propositions together, we have 
proved Theorem 1. 

Appendix E: Computing second cohomology groups 
using generators and relations 

Here we provide some details to justify and explain 
the procedure used for computing second cohomology 
groups in Sec. IV and Appendix G. We focus on the 
f-twisted Z 2n cohomology group H 2 (G', Z 2n ). This in¬ 
cludes H 2 (G',Z 2 ) as a special case (setting n = 1), and 
the treatment for Hf(G', U(l)) proceeds identically, sim¬ 
ply replacing Z 2n by U(l) throughout the discussion. 

The group Hf(G’, Z 2n ) can be computed by finding, 
and distinguishing, all possible equivalence classes of t- 
twisted Z 2n factor sets uj(gi,g 2 ), for gi,g 2 £ G'. Recall 
that such a factor set is any Z 2n -valued function satis¬ 
fying the twisted associativity condition, Eq. (52), and 
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that we are referring to equivalence classes under projec¬ 
tive transformations defined in Eq. ( 54 ). 

Rather than directly studying factor sets, we can 
equivalently study t- twisted Z2„ group extensions of G'. 
Such a group extension is a group E for which Z2 n C E 
is a normal subgroup, satisfying E/^n = G' . An arbi¬ 
trary element e £ E can be written e = au(g), where 
a £ Z2 n , and u(g) is chosen to satisfy 7r[«(<?)] = 3, where 
7r : E —^ G r is the projection map associated with the 
quotient of E by Z 2n . We refer to u(g) as a representa¬ 
tive of g in E. We require the additional property 

u(g)a = a t(9) u(g), (El) 

where t : G' -> Z2 is the twisting homomorphism dis¬ 
cussed in Sec. IIIB. We note that the representative u(g) 
is arbitrary up to projective transformations 

u(g) -> A ~ l (g)u(g), (E 2 ) 

where A (3) £ Z 2n . 

It follows from the definition that 

u(gi)u(g 2 ) = u}{g 1 ,g 2 )u(gig 2 ), (E 3 ) 

where ui(gi,g 2 ) £ Z 2n . Associative multiplication of the 
u(g)’ s, together with Eq. (El), implies that u satisfies 
Eq. ( 52 ), and is thus a f-twisted Z2 n factor set. In addi¬ 
tion, under projective transformations Eq. (E 2 ), the fac¬ 
tor set transforms as in Eq. ( 54 ). So we have shown that 
a group extension is associated with a unique equivalence 
class [w] £ H 2 (G', Z 2n ). 

Now we would also like to show that, given a factor set 
oj (31 ,32), we can construct a corresponding group exten¬ 
sion. We consider a set E whose elements are ordered 
pairs (01,3), where a £ Z 2n and 3 £ G' . We make this set 
into a group by defining the multiplication operation 

(01,31) x (a 2 ,32) = (a 1 a 2 (9l) u;(g 1 , g 2 ), g 1 g 2 ). (E 4 ) 

With this multiplication, it can be checked that E is a 
group, and indeed a f-twisted Z2„ group extension. [ 71 ] 
Choosing u(g) = ( 1 ,3), we have u(gi)u(g 2 ) = 

w(3i, g 2 )u(gig 2 ), as desired. 

It follows from the above discussion that, if we would 
like to construct all possible factor sets (or equivalence 
classes thereof), it is enough to construct all possible 
group extensions. We now describe, in general terms, 
how to do this for a group G' presented in terms of gen¬ 
erators and relations. This procedure is worked out in 
Sec. IV and Appendix G for specific examples. We note 
that in those sections, to simplify the discussion in the 
main text, we slightly abuse terminology and refer to pro¬ 
jective representations, which are group extensions with 
additional vector space structure. This additional struc¬ 
ture is not used in the cohomology group calculations, 
which can be viewed more simply as calculations with 
group extensions. 

To begin, we describe the presentation of G' in terms 
of a finite number of generators hi £ G' ( i = 1 , 2 ,...). 


Note that our goal here is not to define G' abstractly 
in terms of generators and relations, but rather to give 
a description of G' in this manner, assuming that G' is 
already defined by some other means. For every 3 £ G ', 
we choose a fixed canonical form in terms of the genera¬ 
tors, for example 31 = h\h 2 . In general, different choices 
of canonical form are possible for each 3, and fixing the 
canonical form should be viewed as an arbitrary choice. 
Fortunately, while we use the canonical form to justify 
our calculation procedure, it is not necessary to make a 
specific choice in the explicit calculations. It is impor¬ 
tant to note that h is not automatically included as a 
generator, but sometimes it may need to be included, so 
that all 3 £ G' can be written as a product of generators. 

The generators obey a finite number of relations, for 
example, 

h{ = 1 (E5) 

(hih 2 ) 4 = 1, (E6) 

and so on. For the present purposes of general discussion, 
we work in a convention where the right-hand side of each 
relation is the unit element; however, this is not always 
convenient in practice. The relations must be chosen so 
that, given any 31,32 £ G' expressed in canonical form, 
the relations alone can be used to bring the product g±g 2 
to canonical form. 

Now suppose £ is a t-twisted Z 2 n group extension 
of G'. For each 3 £ G' , by making suitable projec¬ 
tive transformations, we can choose a canonical form 
for 11 (g), which is the product of u(hi) corresponding to 
the canonical form of 3 . For example, if 31 = h\h 2 , 
we choose 11 ( 31 ) = u(hi)\u(hz )] 2 , with trivial Z2„ coeffi¬ 
cient. It is always possible to make such a choice, by mak¬ 
ing projective transformations u(g) —> A ~ 1 (g)u(g), where 
A (hi) = 1 . We also choose it(l) = 1 . In addition, if hi and 
h~ x are both generators, we choose u(h~ x ) = [tt(/i,)] _1 , 
which can be accomplished via a projective transforma¬ 
tion A( 3 ) where A( 3 ) = 1 if 3 7 ^ h~ x . 

The relations now become relations for the u(hi), with 
the right-hand side modified to be an arbitrary element 
of Z 2ra , for example, 

[u(h 1 )] 2 = a 1 (E7) 

[u(/ii)u(/i 2 )] 4 = a 2 , (E8) 

for oq, a 2 £ Z 2n . We note that, due to the special choice 
of u(h£ l ) when both hi and hj 1 are generators, we au¬ 
tomatically have OLi = 1 for the relation h t ■ hj 1 = 1 . 
These relations allow us to bring any product u(g\)u(g 2 ) 
into canonical form u(g\g 2 ), up to a Z 2 n phase factor de¬ 
termined by the {oq}. This phase factor is nothing but 
^( 31 , 32 ) S Z 2n , and u(gi)u(g 2 ) = u(g 1; g 2 )u(gig 2 ). The 
set {ai} thus determines u>(gi,g 2 ). We note that the sets 
{cti} can be multiplied according to 

{ai} x {/3 ,} = {aj/3J, (E9) 

which corresponds to the multiplication of factor sets. 
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It is clear that any extension E can be described by a 
corresponding set {ai}. (Note that the converse of this 
statement is not true.) This fact allows us to find all 
equivalence classes [w] £ H 2 {G', Z 2n ) via the following 
procedure. First, we consider the at to be free param¬ 
eters. We then exploit the remaining freedom to make 
projective transformations, where A (g) ^ 1 only if g is a 
generator, to “fix a gauge” for the a,;. After gauge-fixing, 
distinct sets {ai} are inequivalent under projective trans¬ 
formations. Next, we need to determine which sets {ai} 
are consistent, giving rise to an extension E (or, equiv¬ 
alently, to a factor set w). Some sets {ai} can be ruled 
out by algebraic manipulations of the relations; for ex¬ 
ample, one can conjugate various relations by one of the 
generators, which often puts constraints on some of the 
a,. After ruling out some sets {ai} in this manner, one 
can tentatively conclude that the remaining gauge-fixed 
sets {ai} correspond to elements of H 2 {G', Z 2n ). This 
not only gives a computation of the group H 2 (G', Z 2n ), 
but also an explicit parametrization in terms of gauge- 
fixed sets {ai}, with the group multiplication given by 
Eq. (E9). 

To verify this tentative answer for H 2 (G',Z 2 „), one 
needs to show that each {ai} in fact gives rise to a factor 
set. It is enough to do this for sets {ai} that generate 
H 2 (G', Z 2n ). In each case, we can verify the existence of 
the corresponding factor set by, for example, exhibiting 
a projective representation for which the relations realize 
the set {ai}. 


Appendix F: Coupled layer construction 

The fact that all symmetry fractionalization patterns 
are possible on the surface of some d = 3 SPT phase 
(which may be the trivial SPT phase) plays an important 
role in the discussion of this paper. Here, we establish 
this fact using a simple generalization of the coupled layer 
construction of Ref. 37; the discussion in the first part of 
this Appendix closely follows Section IV of that work. 

We consider a symmetry group G , and a fractionaliza¬ 
tion pattern e[u> e ]m[u; m ] (perhaps anomalous) for a d = 2 
state with Z 2 topological order. Here, [w e ]z 2 , [w m ]z 2 £ 
M 2 (G,Z 2 ) are the fractionalization classes of e and m 
particles, respectively. There are no restrictions on [tu e ]z 2 
and [w m ]z 2 ; that is, ([w e ]z 2 , [w m ]z 2 ) is an arbitrary pair 
of elements of H 2 (G, Z 2 ). 

We build a d = 3 system as a stack of d = 2 lay¬ 
ers of SET phases with Z 2 topological order, alternating 
between layers where we label the two bosonic anyons 
as Ei,rrii, and layers where they are labeled e,, M t , as 
shown in Fig. 3. We choose the Ei to have fractional¬ 
ization class [u> e ] z 2 , and the Mi to have fractionalization 
class [w m ]z 2 ■ The e* and rrii have trivial fractionalization 
class. We argue below that such layers can indeed be 
realized strictly in d = 2. 

We assume we have a total of N layers with N 
even, and condense the composite particles Eiei+iEi+% 



FIG. 3. Coupled layer construction. Each layer is a SET 
phase with Z 2 topological order. E and M particles trans¬ 
form non-trivially under the symmetry G, while e and m 
transform trivially. Composite particles indicated by ovals 
are condensed to obtain a d = 3 SPT phase (which may be 
the trivial SPT phase). The particles in dashed boxes remain 
deconfined and uncondensed, and give rise to surface SET 
phases at the top and bottom surfaces. By choosing the frac¬ 
tionalization classes of E and M, surface SET phases with any 
desired symmetry fractionalization pattern can be realized by 
this construction. 


(for i = 1, 3,..., N — 3), and Miim + iMi + 2 (for i = 
2,4,..., N — 2). These particles are bosons with triv¬ 
ial mutual statistics, so they can indeed be condensed 
simultaneously. Moreover, the fractionalization classes 
of these particles are trivial, so they can be condensed 
without breaking symmetry. 

In the state obtained upon condensation, all anyon ex¬ 
citations in the bulk are either confined or condensed. 
Since the symmetry is not broken by the condensation, 
the resulting state is thus a d = 3 SPT phase, which may 
be the trivial SPT phase. At the i = 1 surface, the parti¬ 
cles Ei and mi M 2 remain deconfined, and have fraction¬ 
alization classes [u> e ] z 2 and [w m ]z 2 , respectively. These 
are the quasiparticles of the desired surface SET phase 
with Z 2 topological order and fractionalization pattern 
e[u: e ]iTi[u! m }. The same holds at the i = N surface, for 
the particles En_iCn and Mpj. 

To conclude the discussion, we need to verify that 
the layers in our construction are allowed strictly in 
d = 2. Equivalently, we need to argue that the fraction¬ 
alization pattern e[u e ]mO is non-anomalous for arbitrary 
[w e ]z 2 £ if 2 (G,Z 2 ). To do this, we construct a Z 2 gauge 
theory where the matter field carrying Z 2 gauge charge 
transforms with fractionalization class [w e ]z 2 , and show 
that this gauge theory can arise as a low-energy theory 
for a spin model. 

The Z 2 gauge charge is carried by a multi-component 
boson field b} a , where r labels the sites of a lattice invari¬ 
ant under the symmetry, and a labels the components. 
We take the symmetry operation g £ G to act on the 
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boson field by 


g:bl a ^ r {g) ap b\ rtP . (FI) 


Here, the matrices T(g) are chosen to form a projective 
representation of G whose factor set belongs to the de¬ 
sired fractionalization class [w e ]z 2 . 

We choose a set L of lattice links £ = (r, r') that make 
the lattice into a connected graph respecting the sym¬ 
metry, and introduce a Z 2 gauge field defined on links 
i £ L. On each link £ £ L we introduce a two-dimensional 
Hilbert space, acted on by the Z 2 vector potential er| and 
the Z 2 electric field af. These operators can be thought 
of as 2 x 2 Pauli matrices. Apart from the action of 
space group operations on links, symmetry acts trivially 
on these fields, that is 


9 ■ <y 


The Hamiltonian takes the form 


(F2) 


= (F 3 ) 

ieL p t£p r 


where h,K,u > 0, and the second sum is over a set of 
elementary plaquettes p of the lattice. We may add addi¬ 
tional short-ranged terms consistent with symmetry, but 
will not need to do so for the present discussion. We also 
have to specify the Gauss’ law constraint, which we take 
to be 


= (-!)*■*", (F4) 

r'~r* 

where the product is over those sites r' joined to r by 
some link (r, r') £ L. 

We consider two limits of the Hamiltonian. First, when 
h = 0, the Hamiltonian is exactly solvable, and describes 
a Z 2 gauge theory in its deconfined phase, with gapped 
bosonic matter. The e particles, which are simply the 
bl a bosons, have fractionalization class [w e ]z 2 - To see 
that the m particles have trivial fractionalization class, 
we note that we can integrate out the bosonic matter 
in the limit where u is large, to obtain a pure Z 2 gauge 
theory with gauge constraint n r '~r a 'rr' = 1 - Because 
there is no background Z 2 gauge charge, symmetry acts 
trivially on the m particles, and the m particle fraction¬ 
alization class is trivial. Therefore, this gauge theory 
indeed realizes the e[w e ]m 0 fractionalization pattern. 

We also consider the limit h u, K, which is a con¬ 
fining limit for the Z 2 gauge field. In this limit we may 
put af = 1 , and the gauge constraint becomes 

(_l) b L b ™ = 1. (F5) 

This constrains the number of bosons to be even on each 
lattice site, and defines the Hilbert space for a bosonic 
model, for which the Hilbert space is a product of site 
Hilbert spaces. Because all operators acting within this 
Hilbert space add or remove even numbers of bosons, 


such operators transform linearly under G, which is an 
important requirement for any physical model with G 
symmetry. We thus recover a sensible spin model in the 
confining limit of the gauge theory, and, therefore, the 
gauge theory can arise as a low-energy effective theory of 
such a spin model. We then expect that the deconfined 
phase with e[w e ]m 0 fractionalization pattern can occur 
in this spin model, albeit for some unknown and possibly 
complicated Hamiltonian. 

In certain special cases, it has also been shown via 
construction of exactly solvable spin models ( i.e ., not 
parton gauge theories) that all fractionalization patterns 
e[w e ]?n0 can occur strictly in d = 2. This has been done 
for arbitrary finite, unitary, on-site symmetry [72], and 
also for pAmm square lattice space group symmetry [65] . 


Appendix G: More examples 

1. G = (U(l) x Z 2 ) x pm 


This symmetry is closely related to the case G = U(l) x 
pm , but now with time reversal symmetry added. The 
Zj time reversal forms a semidirect product with U(l). 
The generators are as in Sec. IV C, with the addition of 
the time reversal operation T, and we have the relations 


rp rp rp — 1 rp — 1 

1 x 1 y 1 x 1 y 

= 1 

(Gl) 

TyPxTyiPx 

= 1 

(G2) 

P 2 

X 

= 1 

(G3) 

T X P X T X P X 

= 1 

(G4) 

T 2 

= 1 

(G5) 

TT X 

= t x t 

(G 6 ) 

TTy 

= TyT 

(G7) 

TP X 

= P X T. 

(G 8 ) 


The m symmetry fractionalization is specified by 


rprnrprnrprn—lrpm—1 _ m 

-L x y x y ®txty 

rpm pmrpm— 1 pm _ m 

^ y ^ x ^ y r x u typx 

(P m ) 2 = cr m 

rpm pmrpm pm _ m 

^ X X J- X x ®txpx 

(rrm \ 2 „m 

(I ) = cr T 


T mrpm _ m rpm rpm 

J- x ®Ttx^- x ' 

T mrpm _ m rpm rpm 

1 y ~ <J Tty 1 y ' 

T m pm _ m pmrpn 

*x ^Tpx-^x ' 


(G9) 

(G10) 

(Gil) 

(G12) 

(G13) 

(G14) 

(G15) 

(G16) 


where the a m ’s take values in Z 2 . All the a m, s are invari¬ 
ant under projective transformations of the generators, 
so we tentatively conclude that [w m ]z 2 £ H 2 (G', Z 2 ) = 
(Z 2 ) 8 . To be sure this is correct, we need to show that 
each of the possible 2 8 choices of the a m, s can actually 
be realized by a corresponding factor set. It is enough to 
give a set of projective representations whose cohomology 
classes generate H 2 {G', Z 2 ); this is done in Table II. 
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Rep. number 

T x 

T 

± y 

Px 

U T 

a’s that are —1 

1 

a x 

i 

ia v 

1 

Gpx 5 Gtypx 

2 

cr x 

i 

<J Z 

1 

Gtxpx 

3 

a x 

a z 

1 

1 

G txty 

4 

i 

a x 

<J Z 

1 

Gtypx 

5 

i 

i 

1 

ia y 

<7 T 

6 

a x 

i 

1 

cr z 

artx 

7 

i 

<j x 

1 

a z 

a T ty 

8 

i 

i 

a x 

cr z 

GTpx 


TABLE II. Set of 8 projective representations whose cohomol¬ 
ogy classes generate H 2 (G', Z 2 ) = Z§, where G‘ = pm x Zf\ 
Note that the cohomology classes of the first four represen¬ 
tations listed generate H 2 (pm, Z 2 ) = Z\. The first column 
numbers the representations, 1 through 8. The middle four 
columns specify generators of the group in the corresponding 
representation (time reversal is T = UpK, where K is com¬ 
plex conjugation). All representations in the table are two- 
dimensional. Generators are specified in terms of the Pauli 
matrices a x,y,z . The last column lists those er’s that are equal 
to —1 for the corresponding representation. 


Next, we need to compute H 2 (G', U(l)), noting that 
t(P x ) = — 1 and t(T x ) = t(T y ) = t{T) = 1. Time reversal 
acts trivially on the U(l) coefficients because T is anti¬ 
unitary and T* f2 = f2 2n_1 ; these two effects cancel out 
so that t{T) = 1. We start by specifying 


rptrptrpt —1 rpt —1 _ 

x-L y-*~ x ^ y ™txty 

(G17) 

rpt ptrpt — 1 pt _ 

J -y r x ± y 1 x ~ wtypx 

(G18) 

i P x) 2 = a px 

(G19) 

rpt ptrpt pt _ 

± x r x ± x r x " ^txpx 

(G20) 

X) 2 =a T 

(G21) 

rrl = aptxTir 

(G22) 

tx = ap t xr 

(G23) 

T f P x = appxPX, 

(G24) 


where the a’s take values in U(l). 

Following the analysis of the case of pm symmetry 
without time reversal (Sec. IV C), we adjust the phase 
of Ty to set atypx = 1 (this does not affect apt v ), and we 
can restrict a px ,at xpx £ Z 2 . Next, we can set ap = 1 
by adjusting the phase of T*. Making this adjustment 
modifies ap px -A aTp ap px = a' Tpx , without changing 
other parameters. While this can be absorbed as a re¬ 
definition of ap px , we will keep track of it explicitly, as 
this is important to work out the map p 2 . Next, we 
can conjugate the last three relations by 7", which gives 


UTtx,aTty,u'p px £ Z 2 . Therefore 


rptrptrpt —1 rpt —1 _ 

X-*- y-*~ x y t^txty 

(G25) 

rpt ptrpt —1 pt _ -| 

1 y *x y *x -*■ 

(G26) 

i P x) = a v x £ ^2 

(G27) 

T x P x P x P x = atxpx £ ^2 

(G28) 

X ) 2 = i 

(G29) 

m = [a Ttx G Z 2 ]T*-r 

(G30) 

TX = [«Tt» G Z 2 ]T‘7^ 

(G31) 

m = [a’ Tpx G z 2 ]Pir . 

(G32) 


Note that we have not used the freedom to adjust phases 
of T x or P x . However, adjusting these phases has no 
effect on the a’s. This suggests the result 

H 2 (G',U(1)) = U(1) x (Z 2 ) 5 , (G33) 

with [w]u(!) € H t 2 (G',U(l)) parametrized by [w]u(i) = 
(atxt y ,Oip X ,a tX p X ,apt x ,aTty,a' Tpx ), with the first entry 
a U(l) phase and the last five Z 2 phases. 

To confirm this result, we need to show that each el¬ 
ement is actually realized by some f-twisted U(l) factor 
set. We introduce two-component field operators v r as 
for G = U(l) x pm symmetry in Sec. IV C. We choose 
T x , T y and P x to act on the v r as in Eqs. (74-76), and T 
acts by 


TvrT- 1 = gpvl (G34) 

Here, gp is a 2 x 2 unitary matrix satisfying g^ = 1, 
so that T 2 = 1 acting on v r . We find six families of 
representations, whose factor sets form a generating set 
for H 2 {G\ U(l)) = U(l) x (Z 2 ) 5 : 

1 ■ gtx = gty = g P x = gp = 1 gives a con¬ 
tinuous family of representations with [w]u(i) = 

{atxtyi I ■> I) 1; 1? 1)- 

2. atxty = 1, gty = i, g P x = icr y , gtx = <? z , gr = 1 is a 
representation with [w]u(i) = (1, —1,1,1,1,1). 

8. &txty f; gty f? g px a i gtx @ > gT l is a 
representation with [w]u(i) = (1,1, —1,1,1,1). 

4. a txt y = gty = g P x = 1, gtx = & z , gp = a x is a 
representation with [w]u(i) = (1,1,1, —1,1,1). 

5. atxty = gtx = g P x = 1, gty = v z , gr = cr x is a 
representation with [w]u(i) = (1,1,1,1, —1,1). 

6. atxty = gtx = gty = 1, g P x = v z , gr = cr x is a 
representation with [w]u(i) = (1,1,1,1,1, —1). 

As in Sec. IV C, the above analysis allows us to 
immediately determine the map p 2 : H 2 {G' 1 Z 2 ) -A 
H 2 {G\ U(l)), and we have 


{.atxty■) a px , at xpX ) aptxi apty, ap px ) — 


/ m m m m 
\ G txtyi Vpx') G txpxi G Ttx 5 


m m m \ 
a Ttyi a T G Tpx ) • 


(G35) 
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FIG. 4. Illustration of the operations generating the d = 2 
space group p4mm, the symmetry group of the square lattice. 
T x and T v are translations by one lattice constant along the 
x- and y- axes, respectively. The vertical dashed line is the 
axis for the reflection P x , and the diagonal dashed line is the 
axis for the reflection P xy . 


This implies that, anomaly-negative fractionalization 
patterns are those with a™ = cr™ px = a^ tx = <j™ ty = 
a T a T P x = 1- The group J\f of anomaly-negative vison 
fractionalization classes is thus A f = (Z 2 ) 3 . The disjoint 
sets of SPT phases distinguished by the anomaly test are 
labeled by elements of S = H 2 (G' ,7*2)/N = (^ 2 ) 5 - 


2. G = U(l) x p4mm 

The group p4mm is the space group symmetry of the 
square lattice. We choose generators T x ,T y , T ~ 1 , T ~ 1 , P x 
and P xy . These operations are illustrated in Fig. 4, and 
obey the relations 


t x t v t~ x t~ x = 1 

(G36) 

TyP x T-'P x = 1 

(G37) 

Ty = P xy T x P X y 

(G38) 

Px= 1 

(G39) 

T X P X T X P X = 1 

(G40) 

Pxy = 1 

(G41) 

(PxPxy) 4 = I- 

(G42) 

The m particle symmetry fractionalization is 

by 

specified 

rpmrpmrpm—1 rpm—1 _ m 

x y x y o txty 

(G43) 

rpm pm rpm—1 pm _ m 

^ y r x y r x u typx 

(G44) 

rpm ppm rpm pm 

y -^xy-^x *~xy 

(G45) 

(PD 2 = a™ 

(G46) 

rpm pm rpm pm _ m 

■*- X ^X -*■ X ^X V txpx 

(G47) 

(P m ) 2 = a m 

V xy) pxy 

(G48) 

f pm p>m \ 4 m 

\*x *-xy) ®pxpxy 5 

(G49) 


where the cr m ’s take values in Z 2 . The relation Eq. (G45) 
has no a m parameter, as this can be removed by adjust¬ 


ing the phase T™ —> — T™. The m particle fraction¬ 
alization classes form the group ff 2 (G',Z 2 ) = (Z 2 ) 6 ; a 
generating set of projective representations verifying this 


result is exhibited in Appendix A of Ref. 19. 

To compute the H 2 (G', U(l)) cohomology, 
specifying the relations 

we begin by 

rpt rpt rpt — 1 rpt — 1 _ 

x-L y-*- x -*-y txtxty 

(G50) 

rpt ptrpt— 1 pt _ 

± y r x ± y 1 x ~ u -typx 

(G51) 

rpt _ pt rpt pt 

y ^xy^x^-xy 

(G52) 

( P x ) = a px 

(G53) 

T X P X T X P X = a txpx 

(G54) 

( P X y ) = a pxy 

(G55) 

( P x P xy ) = a pxpxyi 

(G56) 


where the a’s take values in U(l). We note that t(T x ) = 
t(T y ) = 1, while t{P x ) = t(P xy ) = -1. 

First, we adjust the phase of T* to set a typx —> 1. 
In order to leave Eq. (G52) unchanged, we must also 
correspondingly adjust the phase of T x . Next, we adjust 
the phase of P x to set a pxpxy —> 1, which does not affect 
the other relations. Finally, conjugating Eq. (G53) by 
P x , Eq. (G54) by T X P X , and Eq. (G55) by P xy , we have 
ctpx,&txpx,ct P xy G Z 2 . The relations thus take the form 


rpt rpt rpt — 1 rpt — 1 

1 x 1 y 1 x 1 y ~ Oitxty 

(G57) 

rpt ptjit-1 pt _ i 

y x y x 

(G58) 

rpt pt rpt pt 

-L y J -xy- L x^xy 

(G59) 

{Pi? = Upx G Z 2 

(G60) 

T x P x T x Px = a txpx G Z 2 

(G61) 

( Pxy ) = a pxy G Z 2 

(G62) 

(■ PlPly ) 4 = 1- 

(G63) 

This suggests that H 2 (G’, U(l)) = U(l) x 

(Z 2 ) 3 , with 


[w]u(i) G H 2 {G', U(l)) parametrized by [w]u(i) = 

(OCtxtyi &pxj &txpxi Qpxy) • 

To verify this, we proceed as in the case G = U(l) xpm 
in Sec. IV C, and introduce two-component field opera¬ 
tors v r , with r labeling the sites of the square lattice. 
The generators act on the field operators by 

TxVj-T- 1 = ( atxtyY y/2 gtxVr+x (G64) 

P x v r P V 1 = gpxVp^r (G65) 

PxyVrPxy = 9pxyVp^^ rl (G66) 

where a txty G U(l), P x r = P xy r = {y,x), and 

gtx,g px ,gpxy are 2x2 unitary matrices. The action of T y 
follows from Eq. (G59) and is 

T y v r T y 1 = (■ a tx t v y ra:/2 gtyVr+y , (G67) 

where g ty = g V xy9*tx9*pxy 

The following families of projective representations 
form a generating set for H?(G', U(l)): 
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1. gtx = g P x = g P xy = 1 gives a continuous family of 
representations with [w]u(i) = ( cttxty, 1,1, 1). 

2. a t xty = g P x V = 1, gtx = c z , g px = ia y is a projec¬ 
tive representation with [w]u(i) = (1, —1,1,1). 

3- a txty = g pxy = 1, gtx = i& z , 9 P x = o x is a projec¬ 
tive representation with [w]u(i) = (1,1, —1,1). 

4. atxty = gtx = g P x = 1, g pxy = icr v is a projective 
representation with [wWi) = (1,1,1, —1). 

Finally, the map p 2 is given by 

{atxty , ap X , at xpx , a pxy ) — P 2 ([a^m]z 2 ) 

_ / ra m m m \ /p^o\ 

\Ptxty 5 u px 5 °txpx i ^pxy ) • 

Therefore, the anomaly-negative fractionalization pat¬ 
terns are those with cr™ = = = 1- The group 

A f of anomaly-negative vison fractionalization classes is 
J\f = (Z 2 ) 3 . The disjoint sets of SPT phases distin¬ 
guished by the anomaly test are labeled by elements of 
S = H 2 {G',Z 2 )/ U={Z 2 ) 3 . 

We remark that in this case, all the anomalous frac¬ 
tionalization patterns we find can be understood in terms 
of the symmetry U(l) xZf, by choosing different Zf sub¬ 
groups of jAmrn. 


3. G = (U(l) xi Z 2 ) x p4mm 

This is closely related to the case G = U(l) x pAmm, 
but now with time reversal symmetry added. The Z^ 
time reversal forms a semidirect product with U(l). The 
generators are as in Appendix G 2, with the addition of 
the time reversal operation 7”, and we have the relations 


rp T 1 r T ~'—1 rp —1 

1 x 1 y 1 x 1 y 

= 1 

(G69) 

TyP X Ty 1 P x 

= 1 

(G70) 

T y 

— P T P 

— 1 xy- 1 - x- 1 xy 

(G71) 

P 2 

X 

= 1 

(G72) 

T X P X T X P X 

= 1 

(G73) 

P 2 
xy 

= 1 

(G74) 

{PxPxyf 

= 1 

(G75) 

T 2 

= 1 

(G76) 

TT X 

= T X T 

(G77) 

TP X 

= P X T 

(G78) 

TPxy 

a? 

II 

(G79) 


The m particle symmetry fractionalization is specified by 


rpm rpm— 1 rpm — 1 
y x y 

_ m 

®txty 

(G80) 

-rm ppm rpm —1 pm 
y *x y *-x 

_ m 

® typx 

(G81) 

rpm 

1 y 

_ pm rpm pm 

r xy^~ x ^ xy 

(G82) 

{P x f 

_ m 

® px 

(G83) 

rpm pm rpm pm 
-*■ X *~x -*■ X ^ X 

_ m 

® txpx 

(G84) 

{Pxyi 1 

_ m 

Vpxy 

(G85) 

t pm pm\ 4 
\ jr x r xy ) 

_ m 

Apxpxy 

(G86) 

^j~m^ 2 

^.m 
— (J p 

(G87) 

'j-mrj-im 

m rpm'-pm 
~ °Ttx 1 x ' 

(G88) 

rj~m pm 

m pm rpm 
~ °Tpx r x ' 

(G89) 

rpm pm 
’ xy 

_ m pm rpm 

^Tpxy-^xy ' i 

(G90) 


where the a m ’s take values in Z 2 . The m particle frac¬ 
tionalization classes form the group H 2 {G', Z 2 ) = (Z 2 ) 10 ; 
a generating set of projective representations verifying 
this result is exhibited in Appendix A of Ref. 19. 

To compute the H 2 (G', U(l)) cohomology, we begin by 
specifying the relations 


rptrptrpt —1 rpt —1 _ 

-*■ x-*- y-L x ■*-y ™txty 

(G91) 

rpt pt rpt 1 pt _ 

± y r x ± y 1 x ~ LX typx 

(G92) 

rpt pt rpt pt 

^ y r xy^ x r xy 

(G93) 

(P x ) = OL px 

(G94) 

rpt ptrpt pt _ 

± x r X- 1 x r x ~ atxpx 

(G95) 

i^xy) = a pxy 

(G96) 

{P x P xy ) = a pxpxy 

(G97) 

(T 4 ) 2 = ap 

(G98) 

VTl = artxTlY 

(G99) 

T f P x = arpxPir 

(G100) 

P^Pxy = a T P xyP X yY , 

(G101) 

where the a’s take values in U(l). Here, t(T) 

= 1, and t 


is specified for the other generators in Appendix G 2. 

Proceeding first as in Appendix G2, we adjust the 
phase of T y to set a typx —> 1 . In order to leave 
Eq. (G93) unchanged, we must also correspondingly ad¬ 
just the phase of T x . Next, we adjust the phase of P x 
to set a pxpxy —> 1, which does not affect the other rela¬ 
tions. We also adjust the phase of T f to set ap —» 1. 
This modifies ar px —t a' Tpx = a^appx and ar pxy —t 
a' Tp x y = a T X a Tp xy Conjugating Eq. (G94) by P*, 
Eq. (G95) by T X P x , and Eq. (G96) by P xy , we have 
a px ,atx P x,a pX y £ Z 2 . Finally, conjugating the last three 
relations by T* gives a T tx , a ' Tpx , ct Tpxy £ Z 2 . 
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The relations thus take the form 


rptrptrpt —1 rpt —1 _ 

x^- y^- x -L y ^ txty 

(G102) 

rpt pytrpt — 1 T)t _ 1 

L y *x y G x 

(G103) 

rpt p)t rpt pyt 

y ^xy^x^xy 

(G104) 

{P x ) = a px £ ^2 

(G105) 

Px Px T x P x = a txpx £ ^2 

(G106) 

(P X y) = a pxy £ ^2 

(G107) 

(PxPxy ) 1 2 3 4 5 = 1 

(G108) 

(T 4 ) 2 = 1 

(G109) 

rrl = [a Ttx € z 2 ]T 4 T 4 

(G110) 

m = [a’ Tpx G Z 2 ]P 4 T* 

(Gill) 

T^Pxy = [oiTpxy £ ^yV ■ (G112) 

This suggests that Hf(G l 1 U(l)) = U(l) x 

(Z 2 ) 6 7 8 9 10 11 , with 

Mu(i) £ Ht{G’, U(l)) parametrized by 

{.^txtyz t'*p:r ‘ C^txpxi Qpxy: ^ Ttxi ^Tpxi ^‘ > Tp:i:y ) ' 

[ w ]u(i) = 


To verify this, we introduce field operators v r as in 
Appendix G2, for the case of U(l) x p4mm symmetry. 
The action of T x , P Xl P xy and T y is given by Eqs. (G64- 
G67). Time reversal acts by 


lj 1) 1) 1) 1) 1)- 

2- Oitxty 9pxy — 9T — 1> 9tx & i 9px 

ia v is a projective representation with [w]u(i) 

( 1 ,- 1 , 1 , 1 , 1 , 1 , 1 ). 

3- &txty = 9pxy = 9T = 1? 9tx ^ 1 9px 
a x is a projective representation with [w]u(i) 
( 1 , 1 ,- 1 , 1 , 1 , 1 , 1 ). 

d. Ottxty = 9tx — 9px — 9T — f) 9pxy 
ia v is a projective representation with [ctjum 
( 1 , 1 , 1 ,- 1 , 1 , 1 , 1 ). 

5. Oit x ty = g yx = 9pxy = 9tx — @ , gT 

a x is a projective representation with [w]u(i) 

( 1 , 1 , 1 , 1 ,- 1 , 1 , 1 ). 

6. OLtxty — gtx = 9pxy = 1> 9px = @ > gT 

a x is a projective representation with [w]u(i) 

( 1 , 1 , 1 , 1 , 1 ,- 1 , 1 ). 

7 . Oitxty — gtx — 9px — 1 ; 9pxy — @ , gT 

a x is a projective representation with [w]u(i) 

(1,1,1,1,1,1,-1). 

Finally, the map P 2 is given by 


Tv r T 1 =g T v\, 


(G113) 


[ott.xty■) Qpxi Q-txpxi Otpxy t QTtxi Ot-Tpxi ^Tpxy) — P2({p s}) 

_ / _ra m m m m m m m m \ /piiT 

— \ a txtyi G px-> G txpx'> G pxy’> a Ttxi a T G Tpxi G T a Tpxy)’ 


where gr is a 2 x 2 unitary matrix, satisfying g^ = 1 so 
that = 1- 

The following families of projective representations 
form a generating set for U(l)): 

1 -gtx = g P x = g P xy = gT = 1 gives a con¬ 
tinuous family of representations with [w]u(i) = 


Therefore, the anomaly-negative fractionalization pat¬ 
terns are those with cr™ = a™ px = cr™ y = a^ tx = 
a T a Tpx = a T a Tpxy = 1- The group N of anomaly- 
negative vison fractionalization classes is A f = (Z 2 ) 4 . 
The disjoint sets of SPT phases distinguished by 
the anomaly test are labeled by elements of S = 
U 2 (G',Z 2 )/A7= (Z 2 ) 6 . 
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